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ANALYSIS OF GaAs AND Si SOLAR ENERGY HYBRID SYSTEMS 


By 

John H. Heinbockel^ and A.S. Roberts, Jr.^ 

SUMMARY 


Various silicon hybrid systems are modeled and compared with 
a Gallium Arsenide hybrid system. The hybrid systems modeled 
produce electric power and also thermal power which can be 
used for heating or air conditioning. Various performance 
indices are defined and are used to compare the system per- 
formance. The performance indices are: capital cost per 

electric power out; capital cost per total power out; capital 
cost per electric power plus mechanical power; annual cost per 
annual electric energy; and annual cost per annual electric 
energy plus annual mechanical work. These performance indices 
indicate that concentrator hybrid systems can be cost effective 
when compared with present day energy costs. 

Realistic costs and efficiences of GaAs and Si are respect- 
ively $35,000/m^ for 15% efficient solar cells and $1000/m^ for 
10% efficient solar cells. The performance indices show that 
limiting values for annual costs are 10.3 <;:/kwh and 6.8 <:/kwh 
for Si and GaAs respectively. Results demonstrate that for a 
given flow rate there is an optimal operating condition for 
maximiam photovoltaic output associated with concentrator systems. 
Also concentrator hybrid systems produce a distinct cost advan- 
tage over flat plate hybrid systems. 


^ Professor, Department of Mathematical and Computing 

Sciences, Old Dominion University, Norfolk, Virginia 23508. 

2 Professor of Engineering and Associate Dean, School of 

Engineering, Old Dominion University, Norfolk, Virginia 23508. 



LIST OF SYMBOLS 


Si 

= 

silicon 

GaAs 

= 

gallium arsenide 

P 

— 

maximum power 

max 


H3 

= 

efficiency of solar cell 


= 

reference efficiency 

T 

r 

= 

reference temperature 


= 

slope coefficient 

T 

= 

temperature °C 

T 

K 

= 

temperature °K 

E 

= 

solar intensity [w/m^] 

I 

SC 

= 

short circuit current density [amp/cm^] 

V 

oc 

= 

open circuit voltage [Volts] 

FF 

= 

fill factor 

VOCO, ISCO, FIFO 

= 

scale factors 

A. 

1 

= 

constants (i= 1, 2, 3, 4, 5) 

®ELEC 

= 

electric output per square meter of absorber 
area [w/m^] 

^IN 

= 

energy absorbed by system per square meter 
of absorber [w/m^] 

^LOSSES 

— 

loss terms per unit area due to convection 
and radiation losses [w/m^] 

"n 

= 

incident solar flux density [w/m^] 

^THERMAL 

= 

energy removed by coolant [w/m^] 

T 

CK 

= 

temperature (°K) of plate or cell 

'^FK 

= 

temperature (°K) of cooling fluid 

hi 

= 

heat loss coefficient [w/m^ °K] 

"^AIRK 

= 

ambient temperature (°K) 
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e 

= 

emissivity 

a 

= 

Stephan-Boltzmann constant 

• 

m 

= 

mass flow rate [kg/hr] 

C 

P 

= 

specific heat of fluid 

A , 
abs 

= 

area of absorbing surface 

na 

= 

efficiency of flat plate collector 

p . 
1 

= 

i = 1, 2, 3, 4, 5 performance indices 

A 

ap 

= 

aperture area 



A 

c 

TH 

“ 

= theoretical concentration ratio 

^abs 

^WORK 

= 

•5 ( 1 - t^)Qthermal " quality 

\ ^CK / thermal thermal energy 

®LT 

= 

low equality thermal energy 



mC 

K 

” 

— — = effective thermal conductance per 
abs unit area of absorber 

Hi 

= 

efficiency of absorption 

n 

= 

month of year (1 to 12) 

^ELEC 

= 

electric output for nth month per unit area 
of absorber 

^THERMAL^"^ 

= 

thermal output for nth month per unit area 
of absorber 

®WORK 

= 

work output for nth month per unit area 
of absorber 

Ij,(n) 

= 

average solar flux density for nth month 

i 

= 

interest rate 

m 

= 

maintenance cost 

ni 

= 

number of years 

minP4 

= 

minimum value of performance index 4 
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1 . INTRODUCTION 


A typical solar energy hybrid system is illustrated in 
figure 1. The system is a hybrid system in that solar energy 
can be converted to (a) electrical energy by utilizing photo- 
voltaic devices and (b) thermal energy which in turn can be 
used for heating or air conditioning. 

In this report, various types of solar collectors will be 
analyzed so that electric and thermal energy outputs received 
from the collectors can be compared. The photovoltaic devices 
considered for the conversion of solar energy to electric energy 
will be limited to silicon -(Si) and gallium arsenide (GaAs) solar 
cells. The system comparisons will incorporate both energy 
output comparisons and cost comparisons. The cost comparisons 
are based upon current costs and projected costs for Si and 
GaAs solar cells. Comparisons are obtained by defining various 
system performance indices such as capital costs per unit power 
out and annual costs per annual energy out. 

In section 2, mathematical models are constructed to simulate 
the conversion efficiencies of both Si and GaAs photo cells under 
a variety of temperature and light flux conditions. These models 
are in turn utilized in sections 3 and 4 where various collector 
systems are modeled and compared under nominal and annual solar 
insolation simulations. The results of the system comparisons 
are presented graphically and can be found in these later sections. 
Section 3 is limited to system comparisons under nominal solar 
insolations and section 4 deals only with annual solar insolation 
system performance. 

The results of the system performance on an annual basis 
are in turn utilized in section 5 to illustrate various design 
considerations necessary to meet specific power requirements. 

Section 6 discusses the various assumptions that have been 
made throughout the study and section 7 presents a rationale for 
future cost projections of Si and GaAs solar cells. Section 8 
presents the conclusions of the study. 
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The various graphs presented in this report were obtained 
from the computer programs presented in Appendices A and B, by 
utilizing the appropriate conditions in the programs. 

2. MODELS FOR GaAs AND Si EFFICIENCIES 
Linear Model 

Figures 2, 3, and 5 are taken from references 1, 2, and 3, and 
depict the behavior of GaAs and Si solar cells under a wide variety 
of temperatures and light flux intensities. Figure 4 comes from 
H. Hoval and J.M. Woodall, "Optimization of Solar Cells for Air 
Mass Zero Operation and a Study of Solar Cells at High Temperatures, 
which is a quarterly progress report for the period June 1974 to 
October 1974, NASA contract NASl-12812. This report is not a 
readily available reference. If one assumes that the maximum 
power out is directly proportional to the solar intensity E, 
then one can write 


'’max = ^ 


where a(T) is temperature dependent proportionality constant, 
which in general decreases with temperature. From the definition 
of solar cell efficiency 


ns 


efficiency 


. P 

power out _ max 

power in E • (cell area) 


a(T) 

cell area 


one can conclude that the solar cell efficiency depends upon temp- 
erature only. Thus, it is assumed that 

Ti3 = ti 3(T) = n^. [l - 0(T - T^)] (1) 

where is reference efficiency at reference temperature T^ 

and 8 is proportionality constant. 

For Si, B is chosen as .0041, this value of B gives a zero 
efficiency when T - T^ = 243.9° C. For GaAs, B is chosen as .0024 
which gives a zero efficiency when T - T^ = 416.66° C (reference for 



the above coefficients is "On Heat Rejection from Terrestrial Solar 
Cell Arrays with Sunlight Concentration" by L.W. Florschuetz of the 
Mech. Engn. Dept., Arizona State University, Tempe, AZ. This paper 
was received in a private correspondence with the author and is not 
yet readily available) . 


Nonlinear Model 

The assumed solar cell efficiency TI 3 was modeled to conform 
with the experimental data from reference 2 (fig. 3b). The modeled 
efficiency decreases with increasing light intensity. It should be 
noted that cell design will determine the actual behavior of the 
efficiency. In several GaAs cells, efficiency has been observed 
to increase with increasing light intensity (refs. 4, 5, and 6 ). 

For the purposes of this study, the more conservative assumption 
of decreasing efficiency with intensity was utilized in the effort 
to simulate worst characteristics of mass produced cells. 

The following formulas were used to model the efficiency TI 3 
of Si and GaAs solar cells and were derived on the basis of all 
parametric data taken from various sources (refs. 2, 3, 4, 5, 6 , 

7, 8 , 9, 10) . 

I (E,T) • V (E,T) • FF (E,T) (10*+) 

ri 3 = il 3 (E,T) = -55 2 £__ ( 2 ) 

where 1 ^^ = short circuit current density [amp/cm^] , 

= open circuit voltage [Volts] , 

FF = fill factor, 

E = solar flux density [w/m^] , and 
T = temperature (°C) . 

The following empirical relations were used in the nonlinear 
model for efficiency. 
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GaAs 


V„^(E,T) = VOCOg 

+ 7.361(10"3)E *357 _ 2.45(10“3)t 


I (E,T) = ISCO^ 

SC VJ 

( .148E *363 + . 388)T 

(3) 

FF„(E,T) = FIFO„ 

(3 Ca 

+ 1.922E - 1.11(10-3)t 



V^^(E,T) 

FF (E,T) 
s 


VOCOg + |^2.9(10~3)e _ 2.23(10-3)tJ 

ISCOg ^Aq + A^T + A 2 T 2 + AgT^ + A^T't + AsT^^ E/10 

(4) 

FFg(E,T) 


The above equations were derived to fit the data from figures 2 
through 5 and the I^^(E,T) for Si was taken from reference 3. The 
quantities VOCO, ISCO, FIFO are scaling parameters in order that 
the above equations can satisfy the following reference conditions at 
T = 25° C, E = 1000 w/m2: 


GaAs 

V =1.0; FF^ = .85; 

OC G 


I 

SC 


lO^V FF_ 
OC G 


Si 


V 

OC 


.60; 


FF = .85; 
s 


n E 
r 


SC 


10 FF 


OC 


where = efficiency at reference conditions. The constants A^ 

for Si are: 


Ao = .914727 
Ai = .108713(10 -2) 
A 2 = -.695706 (10 
As = .226603(10“^) 
A 4 = .17109(10-5) 
As = -.144039(10 
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Comparison of the linear and nonlinear solar cell efficiency 
models are illustrated in figure 6 for a variety of solar flux 
intensities using a reference efficiency of 15%. 

The models depicted by (1) and (2) will be utilized to represent 
the solar cell efficiency in the following sections. 

3. HYBRID SYSTEM MODELS AND ANALYSIS 

A hybrid system can be visualized in figure 1. A hybrid 
system can produce electrical power and heating or air conditioning. 
Actually a hybrid system produces electric power by photovoltaic 
conversion of the incoming solar flux density and also thermal 
power by absorption of that light flux density which is not 
converted to electrical power. The thermal power can then be 
utilized in heating or air conditioning. Figure 1 illustrates 
a system that could be used for an ammonia-water absorption 
cooling system. 

Various systems for the absorption of the solar flux are 
considered in this report. The various systems are illustrated 
in figure 7 and can be summarized as follows: 

System I : Two flat collectors, one for the collection of 

solar flux to be converted into thermal power and the other 
having silicon (maintained at 25° C) for the production of 
electrical power. 

System II : A single flat plate collector upon which is 

placed silicon solar cells. It is assumed that the temperature 
of the silicon cells is the same as that of the flat plate and 
the silicon efficiency is a function of this temperature. 

System III : A concentrator system having GaAs solar cells 

with efficiency as a function of temperature. The temperature 
can be controlled by a fluid passing through the absorber. 

System IV : A concentrator system having Si solar cells 

with efficiency as a function of temperature. The temperature 
can be controlled by a fluid passing through the absorber. 
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System V ; A concentrator system with no solar cells and a flat 
plate collector with silicon solar cells. The concentrator system 
is for thermal power production and the flat plate collector is 
assumed maintained at 25° C for the production of electrical power. 


Analysis and Assumptions for Systems I and II 

In system I, illustrated in figure 7, it is assumed that 84% 
of the incident solar flux is absorbed by the flat plate collectors. 
The electric power output per unit area of absorber is 

QeLEC = ^3 °IN 


where = -841^^ [w/m^] and na is cell efficiency. 

The energy loss from the collector is 

®LOSS ^ “ "^AIRK^ ~ “ '^AIRK^ 


(5) 


It is assumed that the thermal energy obtained from flat plate 1 
is the energy remaining after losses are accounted for. An energy 
balance produces 


^IN " ^LOSS “ '^THERMAL ® 


or 


“^THERMAL "^AIRK^ '^AIRK^ 


mC 


abs 


— (T - T ) 
' CK FK^ 


( 6 ) 


where 


N 


CK 


AIRK 


incident solar flux (assumed to be 500 w/m^) 
temperature of flat plate [°K] 
ambient temperature [°K] 
emissivity = .04 
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a = Stephan-Boltzmann constant = 5.6697(10“®) w/m^ « 

m = mass flow rate [kg/hr] 

C = specific heat of fluid [kj/kg °C] 

P 

= temperature of fluid entering collector [°K] 

r K 

A , = area of absorbing surface [m^] 


The cost of the flat plate collectors is assumed to be a 
linear function of the thermal efficiency ?i 2 where 


02 


efficiency of collector 


^THERMAL 



and the cost is given by 

cost of collector = 400 ri 2 + 10 

This corresponds to current day costs which range between $53.82/m^ 
($5/ft2) and $645.90/m^ ($60/ft^) for flat plate collectors. 

The silicon solar cell costs are analyzed in section 7 and 
are taken as $1000/m^ for 10% efficient cells and the cost of 
maintaining cell temperatures at 25° C is assumed to have the 
minimal value of $56/m^. This gives a total cost for the flat 
plate 2 of $1056/m^. 

System II is illustrated in figure 7. Assuming that 84% of 
the incident flux energy is absorbed by the solar cells an energy 
balance on system II gives 

^Ij^ “ ^ELEC “ ^LOSSES “ ^THERMAL "" ° 


or 


(1 - n q ) (.84)1 - hi (T - T ) - ea (T** - t“* ) 

N ^ ' CK AIRK-* ' CK AIRK"^ 

mC 

= — (T - T ) 

Aabs CK FK^ 


(7) 
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where 


I 


N 


hi 


T 


CK 


e 


a 


T 

AIRK 

• 

m 

A , 
abs 


incident flux density [w/m^] (500 w/m^ ) 

heat loss coefficient [w/m^ °K] (500 w/m^ °K) 

temperature of cells = temperature of flat plate [°K] 
emissivity (.04) 

Stephan-Boltzmann constant = 5.6697(10“®) w/m^k** 

ambient temperature [°K] 

mass flow rate of fluid 

specific heat of fluid 

area of absorber 


It is known that the second law of thermodynamics gives an 
expression for the maximum thermal efficiency of a heat engine. 


which is the Carnot efficiency, 



where is fluid temperature and is cell temperature. Here 

it is assumed that the fluid will enter collector at a temperature 
Tpj, and exit from collector at a temperature This assumption 

is consistent with the fact that typical large power plants have 
overall efficiencies of 50 to 60% of the Carnot efficiency (ref. 11). 


Actually most solar heat engines go through a Rankine cycle 
which is approximately 50% the efficiency of a Carnot cycle. 
Define 


®W0RK 




^THERMAL 


(8) 


as the high quality thermal energy that can be extracted via 
turbine from the thermal energy received by the solar collector. 

In order to compare the various systems in figure 7, various 
performance indices are defined which will characterize the various 
forms of energy obtainable from a hybrid system. 
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For both systems I and II the following performance indices 
are defined: 


CAPITAL COST 
PEAK (ELEC + MECH) POWER 


(COST FPi 


0 

^ELEC 


)A^bs (COST FP2 


T 


abs 


+ .51 - 


FK 

'CK 


+ COST CELL) A ^ 
abs 


^ ^THERMAL ' '^abs 


CAPITAL COST 
TOTAL POWER OUT 

(COST FPi)A^^^ + (COST FP2 + COST CELL) A 
^^ELEC ^THERMAL^^abs 


CAPITAL COST 
ELEC POWER 

(COST FPi)A^j^^ + (COST FP2 + COST CELL)A^j^g 
®ELEC * ^abs 


where 


A w 
abs 

^ELEC 

^THERMAL 

®WORK 


area of absorber [m^] 

ri3(-84)Ijj 
^^"^CK “ '^FK^ ' 


'^CK / ^thermal 

is the fraction of thermal energy which is 
converted to mechanical work assuming a Rankine 
cycle which is modeled as 50% of a Carnot cycle 
efficiency. 



mC 


K = 




abs 


flow rate parameter 
[kw/m^ °c] 


f 


The same cost figures are used in system II as those presented 
for system I. Systems I and II are summarized in figures 8 and 9. 

Analysis and Assumptions for Systems III and IV 

Systems III and IV are concentrator systems having GaAs and 
Si solar cells respectively and are illustrated in figure 7. 

The following assumptions and notations are used to describe 
and model these systems. 

Let A^p = aperture area, ^^bs ~ absorber area of receiver, 

A 

= theoretical concentration ratio, A = surface area 
abs ^ 

of receiver. It is assumed that — = 2. Then for both systems 

abs 

III and IV the various energy terms per unit area of absorber are: 




where m = optical efficiency of con- 
centrator (assumed to be .7) 


Q 


ELEC 


= risOy where Ha = solar cell efficiency 

J--KT 


■-N 


2h^ ~ '^AIRK^ ~ 2ea (T^^ — 


^LOSSES 


CK AIRK' 


( 10 ) 


mC 

-THERMAL “ A ^ ^'^CK " '^FK^ ~ '^FK^ 

abs 


(11) 


^WORK • ^ \ ^ “ T 


FK 


Q 


CK 


THERMAL 


( 12 ) 


where the notation is the same as that employed to model systems 
I and II. 

An energy balance on the concentrator systems produces 


- Q 




ELEC 


- Q 


LOSSES 


= Q 


THERMAL 


(13) 


13 


I 


and the performance indices i = 1, 2, 3 become 


Pi 


CAPITAL COST 

PEAK (ELEC + MECH) POWER 


(COST CON)A^p + (COST CE;LL)A^^^g 
®ELEC * \bs ®WORK * \bs 


(COST CON) + (CO^ CELL) 
^TH 

^ELEC ^WORK 
C 

TH 


CAPITAL COST 
TOTAL POWER 


(COST CON) + 


(COST CELL) 
C 

TH 


®ELEC ^THERMAL 


'TH 


CAPITAL COST 
ELEC POWER 


(COST CON) + 
^TH 

^ELEC 

C 

TH 


(14) 


The cell costs are taken as $35,000/m^ for 15% efficient GaAs 
solar cells and $1000/m2 for 10% efficient Si cells. Concentrator 
costs are assumed to have a constant value of $156/m2. This high 
cost of the concentrator includes two-dimensional tracking and 
structural supports for wind loads and high concentrations. It 
is a very conservative value when one compares average concen- 
trator costs of $34/m^ from references 12, 13, and 14. Systems 
III and IV are summarized in figure 10. 

Analysis and Assumptions for System V 

System V is illustrated in figure 7 and consists of a con- 
centrator system for thermal energy and a separate flat plate 
system for electrical energy (modeled the same as in system I) . 

Using the same assumptions and notations defined previously, 
the various energy terms associated with system V are: 


14 



(15) 


®IN ^ 1 ^N^TH 


®ELEC na ( • 84) I 


N 


ri3 = solar cell efficiency 


"LOSSES 


2hi + 2ea - TaruK-) 


CK AIRK'' 


(16) 

(17) 


T 


FK 


®WORK \ " "^CK y '‘thermal 


Qn 


(18) 


itiC 


"THERMAL A 


abs 




T ) 
FK^ 


(19) 


It is assumed that the silicon solar cells will operate at maximum 
efficiency and that the energy balance on the concentrator system 
is : 


^THERMAL “ ^LOSSES 

N 


(2o; 


Again, the following performance indices are defined: 


CAPITAL COST 

PEAK (ELEC & MECH) POWER 


(COST CON)A + (COST CELL + COST FP)A 
ap ap 

®ELEC ' ^ap ®WORK ' ^abs 


(COST CON) + (COST CELL + COST FP) 

^WORK 

^ELEC C^jj 


15 


CAPITAL COST 
TOTAL POWER 


(COST CON) + (COST CELL + COST FP) 


Q 


ELEC 


^THERMAL 

*^TH 


CAPITAL COST 
ELEC POWER 


( 21 ) 

concl ' d 


(COST CON) + (COST CELL + COST FP) 
^ELEC 


CONSTANT 


A summary of the equations describing system V are illustrated 
in figure 11. 

Systems I, II, III, IV, and V were compared for a solar 
flux density of I^ = 500 w/m^ at various concentrations and 
flow rates. The computer program used is presented in Appendix A. 

The results of the computer program in Appendix A are 
illustrated in figures 12, 13, and 14. These results utilize 
the nonlinear model for solar cell efficiencies as a function 
of intensity and temperature. Results are for a nominal 
intensity of 500 w/m^. 

In figures 12, 13, and 14 systems I and II have the highest 
capital cost per energy output. All three performance indices 
show capital cost in the neighborhood of $10,000/kw. System V 
has the next highest capital costs — ranging between $3000 to 
$9000 per kw of energy produced. The systems with the lowest 
costs are the concentrator systems with GaAs and Si solar cells 
(systems III and IV) . Capital costs for these systems depend 
upon concentration and operating temperatures and range between 
$900 and $2,000 per kw of energy produced. This is slightly 
above capital costs (installation costs) of nuclear energy 
which is currently running in the neighborhood of $700/kw and 
energy from fossil fuels which costs $550 to $600 per kw 
installed. Gas turbine energy installed costs are around $135/kw. 
However, fuel cost for this type of energy is $2.46 per million 
KJ ($2.60 per million BTU) of energy produced which is expensive. 





Nuclear fuel costs are only 14 per million KJ (IS*? per million BTU) 
of energy produced. The above figures are data obtained from a 
personal communication with C.F. Miller of the Federal Power Commis- 
sion, Washington, DC 20426. In contrast solar energy has no fuel 
costs and depends only upon weather conditions indigenous to the 
area where it is to be utilized. 

Figures 15 and 16 illustrate solar cell electric output as 
a function of concentration for various values of the flow rate 
mC 

parameter K = — — These curves illustrate that at a fixed 

abs 

flow rate the electric power output increases with concentration 
up to a point where the temperature of the cell can no longer 
be maintained at a low value by the cooling fluid. For the 
higher temperatures at the higher concentrations the efficiency 
of the solar cell will begin to rapidly drop off and the electric 
output will go to zero. The maximum power points for the GaAs 
and Si solar cells occur approximately at temperatures of 245° C 
and 160° C respectively for an incident flux density of 500 w/m^ 
assuming a 70% optical efficiency of the concentrator. 

The triangular graphs of figures 17 and 18 illustrate the 
percent energy distribution for the concentrator systems III and 
IV. In these figures 

^LT ~ ^THERMAL “ ^WORK 

where represents the low quality thermal energy remaining 

after fluid has undergone a Rankine cycle to extract useful work 
from the high temperature fluid. 

At each point of the triangular graphs the sum of the ordin- 
ates will add to 100. The ordinates increase from a to A for 
°ELEC' ^ to B for ^ to C for Q^.^. 

At low concentrations there is mostly electric energy and 
low quality thermal energy produced by the concentrator systems. 

As the concentration increases, the temperature rises and the 
useful work that can be obtained from the fluid by a turbine 
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will increase; however, as the temperature increases, the electric 
output is driven to zero. 

Figures 19 and 20 illustrate temperature for various con- 
centrations of systems III and IV. 

4. ANNUAL PERFORMANCE MODELS 


The analysis of systems I, II, III, IV, and V has indicated 
that the concentrator systems III and IV have the lowest capital 
cost. For this reason the following sections will investigate 
the annual energy production of systems III and IV. 

For this purpose, average yearly insolation data from the 
southwest United States has been selected (ref. 15) . These 
average values are illustrated in figure 21 and represent data 
for the simulation of two-dimensional tracking concentrators. 

Let denote the average intensity for the nth month. 

Then the energy balance (13) becomes 


F 


_ (1 - ^'^CK ~ "^AIRK^ 

- 2ca(Tjj^ - T^^^) - K(T^^ - T^^) = 0 


(23) 


If one assumes that the solar efficiency TI 3 is a linear 
function of temperature then (23) can be written as 

^ ^TH' ^FK' ni/ hi, e) = 0 (24) 

For a fixed concentration C one can solve (24) by 

iteration to determine the temperature T . 

CK 

If instead the nonlinear model for efficiency (2) is used, 
then (23) can be written as 
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0 (25) 


F = Ijj(n), ri3(E, Di/ hi, e) = 

where E = n i (n) Again, for a fixed concentration 

one can solve (25) by iterative techniques and determine the 

temperature T . 

CK 

For either the linear or nonlinear model for efficiency, 
ri 3 , corresponding to each month n, we have 

QeleC^’^^ = TliIj,(n)C^jjn 3 (24)DAYM(n) 


^THERMAL " "^FK^ 


^WORK 


~ ^THERMAL 




(24)DAYM(n) 


where DAYM(n) represents the number of days in the nth month. 
Then the annual output from the concentrator system is 


ANNUAL ELEC = 


12 

n?i ^ELEC 


abs 


ANNUAL WORK = 


12 

i5l ^WORK^”^ 


abs 


the capital cost of the concentrator system is given by 


CAPITAL COST = (COST CON)A^p + (COST CELL)A^j^^ 


(26) 

(27) 

(28) 


(29) 

(30) 
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where (COST CON) represents the concentrator costs with two- 
dimensional tracking. The concentrator cost is modeled in two 
different ways. The first representation is 

(COST CON) = CONSTANT = $156/m2 (31) 

and the second model is 

(COST CON) = C^jj + 34 $/m2 (32) 


which is a linear function of concentration. 

The capital costs are amortized over a twenty-year period at 
an 8% interest with an assumed maintenance cost of 2% per year. 
The annual cost can thus be represented as: 


ANNUAL COST 


i 


-ni 

1 - [1 + i] 


+ m 


(CAPITAL COST) 


(33) 


where rii = 20, i=.08, and m = .02. 

Two additional annual performance indices are defined. These 

are 


ANNUAL COST 
ANNUAL ELEC 


ANNUAL COST 
ANNUAL WORK 


(34) 


where ANNUAL ELEC is obtained from (29), ANNUAL WORK obtained 
from (29) and ANNUAL COST is obtained from (33). The equations 
for annual comparison of systems III and IV are summarized in 
figure 21. 

The computer program for the comparison of systems III and IV, 
for annual performance, is given in Appendix B. Graphical results 
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are illustrated in figures 22 through 27. For these figures 
silicon costs were $ 1000 /m^ at 10 % efficiency and galliiim 
arsenide costs were $35,000/m^ at 15% efficiency. 

Figute 22 illustrates performance index P 4 vs. concen- 
tration for various values of the flow rate parameter K. 

Observe that there is a minimum value of ANNUAL COST/ANNUAL ELEC 
for each value of the flow rate parameter. Figure 23 illustrates 
performance index P 4 vs. performance index P 5 for various 
values of the flow rate parameter K. In this figure, note 
that there is a distinct minimum value of ANNUAL COST/ANNUAL ELEC 
for various K values. As the concentration is further increased 
the value of P 5 decreases but under the penalty of increasing 
electric costs. Stated differently, the increased concentration 
raises average cell temperature and reduces electric output. 

The performance index P 5 is not a representative per- 
formance index of true costs as the additional cost of converting 
a high temperature fluid to usable work has not been added to 
the capital costs. Approximate additional cost for converting 
this energy is $333/kw for a high temperature fluid. [One 
possible model for this additional cost would be 333/ - T„^) 

where T is temperature of cells and is fluid tempera- 

CK r K 

ture leaving heat engine.] 

Figures 24 and 25 depict annual energy outputs from the 
concentrator systems III and IV. Figures 26 and 27 illustrate 
solar cell cost vs. minimum value of P 4 for various flow 
rate parameter values of K. Listed alongside these curves 
are approximate concentrations where minimum values are obtained. 

5. DESIGN CONSIDERATIONS 

A hybrid system used in conjunction with conventional heating 
and air conditioning equipment must be designed to carry a cer- 
tain percentage of the daily load requirement. The following 
are calculations to give a first estimate for cost and sizing 
of a hybrid system to supply energy to a residential house. 


The following assumptions will be made: 

1. Heating requirements: = 69.9 ( 10 kwh/yr 

= (25 BTU/hr ft^) (265 days/year) (1500 ft^) 

2. Air conditioning requirements: = 26.4(103) kwh/yr 

= (25 BTU/hr ft^) (100 days/year) (1500 ft^) 

3. Electrical power demand: = 8(103) kwh/yr 

4. Area of house: = 139 m^ = (1500 ft^) 

5. Cost of GaAs: = $35000/m^ 

At a concentration ratio of 400, one can assume the following 
energy values per square meter of aperture area (values taken 
from fig. 25) 

ANNUAL QgLEC “ kwh/m^ yr 

ANNUAL QtherMAL " kwh/m2 yr 

assume total thermal load is 69.9 (10^) kwh/yr and total electric 
load is (26.4 + 8)(103) = 34.4(103) kwh/yr. 

Let A denote aperture area and A , denote absorber area 
ap ^ ab 

(photocell area) . Then to meet the thermal demand one would 
require 


1700 kwh/m2 yr • A^^ = 69. 9 (10 3) kwh/yr 


or 


A = 41.12 m2 = (442 ft2) 
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is the aperture area required to meet this demand. This aperture 
area would supply the following electrical power 

= 300 kwh/m^ yr • 41.12 m^ = 12.34 kwh/yr 

which is 35.8% of electrical energy requirements. The required 
solar cell (absorber) area necessary is 

A 

^ V. = = .1028 m2 = 1.106 ft2 

ab 400 

The approximate cost of such a system can be divided as in 
table 1. 

The cost of such a system depends upon concentrator costs 
which could be major costs of the system. There is currently 
no reliable data on concentrator costs and if one takes the 
cost figures as representing lower and upper bounds for system 
cost, an average cost would be $8738.47. 

Of course this is only a rough estimate of system cost. There 
will be economic variations with respect to geographic location and 
weather conditions. It is envisioned that such systems will be in 
widespread use in commercial and industrial applications within the 
next 25 years. 

6. DISCUSSION OF MODELS, ASSUMPTIONS, AND RESULTS 


The models constructed in this report are for steady-state 
operating conditions of the systems under investigation. The 
models represent average operating conditions of the various 
systems and do not consider transient responses. 

The concentrating devices are assumed to give uniform solar 

mC 

cell illvimination and K = — represents an effective uniform 

abs 
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thermal conductance per unit area of absorber. Values of K between 
.01 and .1 represent passive cooling such as wind cooling. In com- 
parison, values of K between .1 and 10 represent active cooling 
with high mass flow rates. In the report K is the heat removed 
from the absorber area per unit of absorber (cell) area per unit 
temperature difference between coolant inlet and outlet temperatures. 
Via a secondary loop this heat is available to drive a vapor engine 
at an assumed 50% Carnot efficiency operating between some tempera- 
ture limits. Heat exchanger losses are ignored. 

There is a wide spectrum of costs in the solar field for flat 
plate and concentrator collectors. There is also a wide spectrum of 
costs for solar cell devices utilizing silicon and gallium arsenide. 
The cost values used in this study were chosen as conservative values 
under the assumption that the costs of $1000/m2 for 10% efficient Si 
and $ 35000 /m^ for 15% efficient GaAs solar cells together with concen- 
trator costs are greater than "other" siibsystem costs. This assump- 
tion is valid for concentrator costs between $30/m2 ^^d $150/m^ and 
concentrations ranging between 230 and 1200 for GaAs systems and 
concentrations between 6 and 30 for Si systems. 

There is not a great deal of information in the literature on 
concentrator costs for full 2-D sun tracking concentrator devices 
and so throughout most of this study the concentrator costs were 
assumed constant ($156/m^) (refs. 12, 13, 14, 16, 17, 18, 22, 23, 

24, 25). In figure 26, the concentrator cost was assumed to vary 
linearly with concentration (34.0 + 0.122 Cn^„)$/m^. 

Ill 

The annual energy produced by the concentrator systems was 

divided into electrical energy, high quality thermal energy, and 

low quality heat. The high quality thermal energy was that energy 

that can be extracted from a high temperature fluid via a vapor 

engine with an assiamed 50% Carnot efficiency. The annual values 

of electrical energy and high quality thermal energy (Q.^^^^) are 

VVUKJx 

illustrated in figure 25. The cost models do not consider the 
additional cost of converting the high quality energy to useful 
work. A more detailed cost analysis would be necessary for such 
a comparison. 
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There are some differences in the literature as to how solar 
cells behave under high concentrations of solar flux energy. In 
this report it was assumed that solar cell efficiences decreased 
with increased illumination. Some investigators report that the 
opposite is true — that is, solar cell efficiencies increase with 
increased illijmination (refs. 4, 5, 6). The solar cell design 
will determine the actual behavior of efficiency. If cell effi- 
ciencies increase with increased illumination, then the results of 
this study can be taken to represent very conservative estimates 
as to how the various systems perform and the cost estimates given 
in the figures of results must be lowered. Another area where no 
data is available is lifetime performance of solar cells under high 
flux densities. 

Concentrators with high concentration ratios 100 to 10,000 
have been developed for the space program (ref. 19) . The technology 
exists for constructing high concentration systems but little data 
is available on costs for such concentrators. 

7. COST PROJECTIONS FOR GaAs SOLAR CELLS FOR 
TERRESTRIAL APPLICATIONS 


The hybrid system performance and economic analysis has been 
undertaken with uncertain component cost information, but no apolo- 
gies are necessary. The intent was to find the limits of performance 
and to establish comparisons between GaAs and Si solar cell hybrid 
systems. It was realized from the outset that solar cell costs — 
especially for Ga^_^Al^As-GaAs cells — would be pivotal, and system 
performance results are reported over parameter domains where this 
assumption remains valid. The cells of interest, grown by the 
liquid-phase epitaxy process, are currently "hand-made" for experi- 
mental purposes; though requiring little material per wafer, their 
cost is understandably high because of the labor (skilled labor) 
intensiveness of the processing. Cost reductions must ultimately 
be predicated on market development for photovoltaic devices and 
mass production techniques. 


Even while costs remain high for the GaAs hetero junction cells 
there are compelling advantages relative to Si cells which motivate 
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continued basic research and stimulate system feasibility studies. 
Among other researchers, Stevenson in the proposal report, D.A. 
Stevenson, "Thin Film Gallium Arsenide for Low Cost Photovoltaic 
Solar Energy Conversion," Report No. CMR-P-73-17, Center for Mater- 
ials Research, Stanford University, Stanford, CA, 1973 (not readily 
available), has pointed out these advantages: 

1. The bandgap is a better match to the solar spectrum, 
therefore better efficiencies can be obtained. 

2. The bandgap is direct, therefore considerably thinner 
cells can be used to absorb the solar energy. 

3. The bandgap is greater, therefore (a) the power output 
with increasing temperature is greater, and (b) the output vol- 
tage is greater. 

4. The minority carrier lifetimes and diffusion lengths 
are less . . . . 

Stevenson also argues convincingly concerning the natural 
abundances of the prime materials, gallium and arsenic. Materi- 
als are available for large quantities of thin (5 to 100 ym) 

GaAs cells; there is needed only the stimulus to reduce the 
labor-intensive current manufacturing processes. 

It will be the intent of the remaining portion of this 
section to describe the current cost basis used in this report 
for Si cells and for' GaAs hetero junction cells. Also a scenario 
is outlined depicting a plausible cost reduction projection for 
GaAs cells based on an expanding market for photovoltaic electric 
power generation. 

In 1975, a few materials laboratories would quote a price 
for small quantities of experimental hetero junction GaAs cells, 
and the number was in the neighborhood of $800,000/m^, a staggering 
figure to face for the photovoltaic/ thermal power system designer. 
Prominent researchers in the field were pragmatically uninterested 
in predicting where current costs might go should a market develop. 
A number of people in the solid-state electronic device industry 
were asked if they would be willing to extrapolate their experi- 
ence with LED and solid-state laser developments over to GaAs 
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solar cell production, and estimate a de-escalating cost curve 
with growing demand — most were not willing to do this, including 
H. Kressel of RCA Laboratories and M.B. Panish of Bell Laboratories 
(ref. 30). The reason given was basic dissimilarity of the devices. 
However, it was learned that two different producers of LED devices 
were employing epitaxial structures. They were currently 
growing an epitaxial layer on a GaAs substrate, dicing-up the 
wafers, and providing electrical contacts at a cost of about 
$3.50/cm^ ($23/in.^). This yields a "current" cost estimate for 
GaAs cells equivalent to $35,000/m^ assuming the technology extrap- 
olates to solar cell production in large numbers. Although it 
has no impact on "current" cell costs, the continuous cell growth 
processes which are being studied (ref. 26) , must inevitably con- 
tribute to diminishing GaAs solar cell costs. 

Solar cell array manufacturing costs (exclusive of substrate 
and encapsulation costs) for "current" Si cells for terrestrial 
application are variously quoted in recent literature: $2000/m^ 

(refs. 27 and 15); $1250/m2 (ref. 28). For the sake of this study, 
however, a value of $1000/m2 is used (personal communication with 
Gilbert H. Walker of NASA Langley Research Center, Hampton, VA 
23665) , based on the knowledge that Solarex Corporation of Rockville, 
MD is producing commercially small Si cell arrays at a cost of 
$10/peak watt. For a 10% efficient cell this is equivalent to the 
$1000/m2 cost value which is used in the analysis. 

A scenario is outlined here demonstrating how GaAs cell costs 
might ultimately reach a competitive level; the scenario is based 
on the work of Wolf (ref. 28) and Baum (ref. 29) combining an elec- 
tric power market penetration projection with the past cost history 
of Si solar cells. Wolf sees the rising costs of conventional 
sources of power and growing public awareness as setting the stage 
for rapid power generation market penetration by photovoltaic sys- 
tems. Significant prototype application is envisioned to commence 
around 1980 providing electric (and thermal) power for buildings, 
residential and commercial. From figure 3 in the paper by Wolf 
(ref. 28) , table 2 is constructed showing anticipated annual 
electric energy from photovoltaic units on buildings. These 
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values are converted to average kilowatts for a year and then to 
"peak" kilowatts or installed capacity by assuming a plant capacity 
factor of 50%. 

Wolf has argued the plausibility of this level of penetration 
over the next 50 to 60 years. Although the numbers appear large 
they represent a small percentage of projected United States total 
energy consumption in the given years. For later reference atten- 
tion is called to the value for installed or "peak" kilowatts for 
the year 1990, a point in time representing the starting date for 
rapid commercial addition of photovoltaic equipment. 

Although GaAs solar cell fabricators have been reluctant to 
predict future cell costs, VARIAN of California has released cell 
fabrication "experience" factors which gauge costs of epitaxially 
grown cells as a function of peak photovoltaic kilowatts installed; 
this is based on their work with GaAs cells under sunlight concen- 
tration (unpublished data from R.L. Bell VARIAN Co., Palo Alto, CA) . 
Based on the VARIAN projections a market growth to 10® kw peak 
installed will be required to achieve a 100-fold decrease in current 
GaAs cell costs. This kilowatt level coincides with the year 1990 
according to the projections of Wolf. 

Historically the unit cost of production falls off exponen- 
tially with rapid market penetration (refs. 27 and 29). If the 
Si solar cell costs of Baum (ref. 29) are plotted along with the 
"current" cost figure of $10/watt (the Solarex Corporation cost 
derived assuming a 10% cell and a one-sun value of 1 kw/m^) , a 
semi- logarithmic slope of -.23 yr is found. This curve is shown 
in figure 28. For GaAs hetero junction cells the "current" cost is 
$35,000/m2. If the number is modified assuming a 15% efficient 
cell and a one-sun value of 1 kw/m^, a dimensional conversion yields 
a current cost of $233/watt, compatible with the data for Si solar 
cells. Extrapolating from this cost using the slope of -.23 yr , 
a value is derived in figure 28 for the point in time where the 
GaAs cell cost will have dropped by a factor of 100 as suggested 
by VARIAN. Fortuitously, this date is 1995, coinciding approxi- 
mately with the date of significant market penetration (1990) 
suggested by Wolf and being consistent with the level of installed 
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kilowatts necessary to stimulate cell production as indicated 
by VARIAN. 

In terminating this section of the report the following 
conclusions are offered: 

1. With significant prototype testing commencing by 1980, 
installed (peak) photovoltaic capacity could reach 10® kw by 
1990 (Wolf, ref. 28). 

2. A 100-fold cost reduction for GaAs cells is feasible as 
the market achieves 10® kw installed (VARIAN) . 

3. If silicon cell cost reductions observed from 1958 to 1975 
are used as a gauge and the same semi-log slope (-.23 yr “^ ) is applied 
to "current" GaAs cell costs, a 100-fold reduction is predicted to 
occur by 1995 concurrent with Wolf's market projections to 10® kw 
installed . 

4. Convergence of Si and GaAs cell costs are indeterminate 
at this time looking ahead 20 to 40 years. Convergence may be 
possible, however, because GaAs cells (a) require less material 
per unit cell area and (b) can display higher efficiencies rela- 
tive to Si cells, especially under the desirable condition of 
sunlight concentration. 


8 . CONCLUSIONS 

Current costs and efficiencies for GaAs and Si solar cells are: 
$35,000/m^ for 15% efficient GaAs cells and $1000/m^ for 10% effi- 
cient Si cells, both efficiencies at AMI. 

Limiting values for annual energy costs from GaAs and Si con- 
centrator systems have the following range of values : a GaAs 

concentrator system ranges between 2<? and 6.8<: per kwh and a Si 
concentrator system ranges between 2.5<= and ll*? per kwh. The 
ranges in annual energy costs reflect the different assumptions 
on concentrator costs which include full two-dimensional tracking. 

For a given flow rate, there is an optimal operating condition 
for maximum photovoltaic output of both GaAs and Si hybrid systems. 
This can be seen by examining figures 15 and 16. 
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The high concentration hybrid systems offer a distinct cost 
advantage over flat plate hybrid systems because the concentration 
increases solar flux density and decreases the solar cell area. 

Solar hybrid systems for the heating and cooling of buildings 
seem to be economically competitive with existing energy sources. 
Additional cost studies should be undertaken to calculate "total 
system costs" as this study did not figure in costs of cooling 
equipment, hot water storage for heating, or turbine energy con- 
version costs. 

As the cost of solar cells decreases, optimum system perform- 
ance from Si and GaAs hybrid systems can be achieved at lower con- 
centrations. For Si cell costs of $50/m^ (NSF goal) and GaAs a 
factor of 20 more expensive, optimal Si performance can be achieved 
at concentrations around 10 while optimal GaAs performance can be 
achieved at concentrations of around 100. 
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Table 1. Approximate cost for solar photovoltaic system. 


1. GaAs solar cell cost: 

(35000) (.1028) 

2. Concentrator cost with tracking: 

( $34/m^) (41. 12 m^) lower bound 
or 

($156/m2) (41.12 m^) upper bound 

3. Battery storage cost: 

(COST/kwh) (1/EFFIC) (HRS OPERATION) (A^j^) (PEAK FLUX OUTPUT) = 
($40/kwh) (1/.65) (3 HRS) (.1028 m2) (8 kw/m2) 

4. Conversion cost: 

(COST kw) (PEAK FLUX OUTPUT) (A^j^) = 

($100/kw) (8 kw/m2) (.1028 m2) 

5. Hot water storage costs: 

($. 2642/litre) (3785 litres) = ($l/gal) (1000 gal) 


= $ 3598.00 

= 1398.08 

= (6414.72) 

= 151.83 

= 82.24 

= 1000.00 


TOTAL COST 


$ 6230.15 

($11,246.79) 


Table 2. Solar /photovoltaic market penetration 


Year 

Electric Kilowatt-hours 
per year 

Average 

Kilowatts 

Installed 

Kilowatts 

1990 

5.0 

X 

109 

5.7 

X 

105 

1. 1 

X 106 

2000 

2.6 

X 

10^1 

3.0 

X 

10^ 

6.0 

o 
1— 1 

2010 

1.0 

X 

roi2 

1.1 

X 

108 

2.2 

X 108 

2020 

1.8 

X 

1012 

2.0 

X 

108 

4.0 

CO 

o 

1 — 1 

X 

2030 

2.8 

X 

1012 

3.2 

X 

108 

6.4 

X 108 

2040 

3.5 

X 

1012 

4.0 

X 

108 

8.0 

X 1Q8 
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Figure !• A solar energy hybrid system 
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Fig. 3. Relationship between for silicon 

PB and temperature for E = var (same 
symbols as in Fig. 2) . 
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Fig. 5. Efficiency of silicon PB as function 
of temperature for E = var (designa- 
tions the same as in Fig. 2) . 


Figure 3. Si solar cell characteristics from reference 2 






EFFICIENCY 


TEMPERATURE, °C 

Figure 3. Ga, A1 As-GaAs solar cell parameters as a 

] "X X 

function of temperature using a Zenon light 
source as a solar simulator. Device effici- 
encies were calibrated with the solar simu- 
lator at JPL. 


Figure 4 . GaAs characteristics from Progress Report 
cited on page 5. 
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Fig. 3. Short-circuit current output as a function of Fig. 4. Open-circuit voltage output as a function of 
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Figure 5a. Si performance from reference 3 
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Figure 6 


Comparison of models for Si and GaAs 
efficiency. 




SYSTEM I 


SYSTEM II 


FLAT PLATE COLLECTOR (THERMAL FLAT PLATE COLLECTOR 

POWER) SEPARATE FLAT PLATE WITH SILICON 

COLLECTOR (Si) (ELECTRIC .POWER) 



SYSTEM III 


SYSTEM IV 


CONCENTRATOR WITH GaAs 


CONCENTRATOR WITH Si 



CONCENTRATOR FOR THERMAL POWER FLAT PLATE 
COLLECTOR WITH Si FOR ELECTRIC POWER 


Figure 7. Solar collectors for hybrid systems. 



SYSTEM I 


FLAT PLATE COLLECTOR FOR THERMAL ENERGY 
FLAT PLATE COLLECTOR (SILICON SOLAR CELLS) FOR ELECTRICAL ENERGY 



ENERGY BALANCE (SILICON PLATE) 


^I^ ■ ^ELEC “ ^LOSSES " ° 

(1 - ri3)(.84)I^ - hi(T^^ - T^Irk) - "^(^CK ' ^yllRK) = ° 


ENERGY BALANCE FOR THERMAL ENERGY 

^I„ ■ ^LOSSES ■ ^THERMAL " ° 

N 

4 4 

^THERMAL " ’^^^N ' ~ ^AIRK^ " " ^AIRK^ "" A^^g " ^FK^ 

COST CON = 400ti 2 + TO $/m2 na = EFFICIENCY OF COLLECTOR = ^^HERMAL 


COST FLAT PLATE PLUS SILICON - $1000. + 56. 


Figure 8. 


Summary of equations describing system I. 



SYSTEM II 


FLAT PLATE COLLECTOR FOR THERMAL ENERGY HAVING SILICON 
WITH EFFICIENCY OF SILICON AS FCN OF TEMPERATURE 



ENERGY BALANCE 


^I,, ■ ^ELEC " ^LOSSES " ^THERMAL ° 
N 


mC 


(1 ~ n3)(«84)I|^ - hi(Tp|^ - T^jpi^) - ea(Tji^, - T^jpi^) = ^ 

PERFORMANCE INDEX P 


N ■■^'■CK " 'AIRK^ " ^^^'CK “ 'AIRK' “ A^ ^^CK ‘ ^FK^ 


Pi = 


CAPITAL COST _ 


PEAK(ELEC + MECH) POWER 


^ELEC ^abs ^ ^THERMAL ’ ^abs 


CAPITAL COST _ (^OST FP, + COST CELL)A^^^^ 


^2 “ total power out 


^^ELEC ^ ^THERMAL^^abs 


P, = 


CAPITAL COST ^^OST FPjAabs + (COST FP^ + COST CELDA^^^ 


3 ELEC POWER 


^ELEC \bs 


Figure 9. Summary of equations describing system II. 
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III - GaAs SOLAR CELLS 

SYSTEMS III AND IV 

IV - SILICON SOLAR CELLS 


^LOSSES 



COST CON = 100 + ni80. $/m2 

[INCLUDES TRACKING STRUCTURAL 
(WIND) SUPPORTS] 

ENERGY BALANCE 

^i|^ ■ ^ELEC “ Bosses ^ *^thermal 

PERFORMANCE INDEX P 


A,„ = APERTURE AREA 
ap 

A . , * ABSORBER AREA OF RECEIVER 
abs 

Ag = SURFACE AREA OF RECEIVER 
A 

C,u = = THEORETICAL CONCENTRATION RATIO 

"abs 

'^abs 

'^If^ " ’^l^N^TH ^ELEC " 

^LOSSES " ■ ^AIRK^ " ^AIRK^ 

mC 

^THERMAL " A^ ^^CK " "^FK^ 


Pi = 


CAPITAL COST 
PEAK(ELEC + MECH)POWER 


(COST CON)A^ + (COST CEll){A^^^) 


TLEC 


A. . + 

abs 


ap 

.5(1 - ; ek ) 

V 'CK^ 


^THERMAL 


abs 


COST CON + ^OST CELL 


^ELEC " T^) ^THERMAL 


COST CELL 

_ CAPITAL COST _ ^TH 

2 - TOTAL POWER ' Q^leC " ^THERMAL 


COST CON + ^OST CELL 
n _ CAPITAL COST _ ^TH 

^^3 “ ELEC POWER 


Figure 10. Suiranary of equations describing systems III and IV 
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SYSTEM V 


CONCENTRATOR FOR THERMAL ENERGY AND SEPARATE FLAT PLATE 
COLLECTOR FOR ELECTRIC ENERGY WITH SILICON SOLAR CELLS 



ENERGY BALANCE 

^THERMAL " ^Ij^ " ^LOSSES 


COST CON = $100 + ni80. 

COST FLAT PLATE + SILICON = $1056/m2 


PERFORMANCE INDEX P 

CAPITAL COST 


1 


PEAK(ELEC + MECH) POWER 


(COST CON )A^p + ( COST CELL + COST FP)Agp 


%EC * ^ap " T 


FK 


CK> 


^THERMAL * abs 


COST CON + COST CELL + COST FP 
n . r(i ^FK\ ^THERMAL 


P9 = 


CAPITAL COST COST CON + COST CELL + COST FP 


TOTAL POWER 


n . ^THERMAL 
^ELEC C^j^ 


n _ CAPITAL COST _ COST CON + COST CELL + COST FP _ 

P 3 - • E T e ' C ' p o wer CONSTANT 


Figure 11. Summary of equations describing 
system V. 
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SYSTEM IV 


SYSTEM III 


NUCLEAR 


K 

FOSSIL 
1975 COSTS 


GAS 

TURBINE 






CAPITAL COST 
ELEC POWER 


1 X 105 



It 

or> 



CONCENTRATION 


Figure 14. P 3 vs. concentration 
















ANNUAL ENERGY PRODUCTION 


>- 
. <c 
3: o 

3 I 


8.8 


9.3 


9.8 9.8 9.7 Q r 


7.8 


| ?-0 , 8.8 


8.0 


7.3 7.2 


JF MAMJJ ASOND 


I,^,(n) = AVG. INTENSITY FOR nth MONTH 


(1 - “ ^airk^ " ^ 


‘CK 'AIRK 


CK ‘AIRK' 


'CK ‘FK' 


Calculate 


^ELEC^"^ " 

nilN(n) 

C^^H 3 ( 24 )DAVM(ri) 

^COOL^”^ " 

- 

Tfk> 

^WORK^"^ " 

.5 ^1 ■ 

■ T^) ‘“cool'"’ 


/12 

\ 

ANNUAL ELEC 

’ 

\n=l 

‘JeLEC<"> ) "abs 


12 


ANNUAL WORK =( S ^WORK^"^)^ 


. n=l 


abs 


( 


CAPITAL\ /C0ST\ /„ ^ . /COST \ 

COST J "UON / ^ ap^ \CELL / 


Wabs) 


ANNUAL COST = 


(1 <■ D" • i 
(1 + i)" - 1 



/CAPITAL^ 
V COST / 


Performance Index 4 = P = COST Performance Index 5 = P = COST 

Performance index 4 - P^ ANNUAL ELEC Performance index b - P 5 ANNUAL WORK 

i = .08; n = 20 years; m = maintenance cost = .02 


Figure 21. Stunmary of equations for annual 
output from systems III, IV. 
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ANNUAL COST 
ANNUAL ELEC 


I 



CONCENTRATION 

Figure 22. Performance index P 4 vs. concentration 
for various flow rate parameters K. 
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ANNUAL ELEC 







0 


con 


Figure 25. Annual energ 
of aperture 


ANNUAL ELECTRIC PLUS THERMAL ENERGY 




K = 3 


K - .75 GaAs 



ANNUAL 

ELECTRIC 

ENERGY 





SOLAR CELL COST [$/m2] 


100,000 


\450 

475ii) 


210 350 



CONCENTRATION 

e 325"^ 
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NONLINEAR MODEL 


LINEAR MODEL 
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Figure 26. Solar cell cost vs. min P 4 for concen- 
trator costs a linear function of 
concentration . 
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APPENDIX A 

COMPUTER PROGRAM FOR COMPARISON OF 
SYSTEMS I, II, III, IV, V 
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€FFIC [i»3GXB,F0R 


FORTRAN V.4(21e() /KI 


15-.SEP-75 


9:05 


PACK 1 


B030I 


SUBROUTINE EFFICfVOCOG, ISCOG. FI FOG. VOCOS, 1 SCOS.fi FOSJ 

00302 


DIMENSION A(6) 

30303 


READ ISCOG. ISCOS.ISCG.ISCS 

00304 


NRITEI3. 100) 

00305 

103 

FORMATf ■ 1 ' ) 

. 30036 


Atl)=500. 

30007 


A(2)=1000. 

3030B 


.A1.U=500_0. . . 

30309 


AC4)=10000. 

•30310 


A(5)=53000. 

30311 


A(6)=100000. 

00312 


00 1 1=1,6 

30313 


E=A(I) 

_^3 031.4— 


WR1TE(3.101) E 

00015 

101 

FORMAT!//, 1 X , ' INTENSITY = ' ,F9.1 ,//T3. 'TEMP' Til 'EFFG'T21 , ' VOCOG 

00316 


2131,' ISCG'T4I, 'F1DLG'T51 . ' EFFS ' T61 . ' VOCOS ' T7 1 , ' ISCS'TPl . 

00017 


3'FILDS' .//) 

00318 


DO 2 J=1 ,13 

00019 


TEMP=(J-1 )»25. . 

33320 


CALL CH AR (E, TEMP, TSCG.VOCG, FI LUG, FFFG.l .VOCOG, TSCOG.FIFOGl 

00321 


IFCTEMP .GT. 200„) GO TO 3 

30322 


CAUL CriARlE.TEMP, tSCS.VOCS, FIDOS, EFFS. 2, VOCOS, ISCOS. FIFOS) 

30323 

3 

CONTINUE 

30324 


' IFCTEMP .Gi. 200.)“ GOTO 4 

00025 


MR ITS (3,12) TEMP.EFFG, VOCG,ISCG,FILt.G,EFFS.VOCS.ISCS,FlU,S 

. 03326 

12 

10RMAr(9(lX,F9.5)l 

00327 


GO TO 2 

00328 


4 NRITE(3.12) TEMP , EFFG , VOCG, ISCG, F IDI.G 


30329' 2 

CONTINUE 


' 00030 1 

CONTINUE 


00331 

RETURN 


30332 

END 



SUBPROGRAMS CAOOED 
CHAR 


SCALARS 

AND 

ARRATS [ "♦ 

" NO 

EXPLICIT DEFINITION 

- "%■’ 

NOT 

REFERENCED 1 




♦VOCOS 

1 

ISCOG 

2 

♦ VOCS 

3 

*KFFG 

4 

*E 

5 

4V0C0G 

h 

♦ EFFS 

7 

♦FIFOS 

10 

♦FILLS 

11 

♦j 

12 

A 

I 3 

.S0001 

21 

.S0003 

22 

ISCG 

23 

♦FIFOG 

24 

ISCOS 

25 

♦FILLG 

26 

*TEMP 

27 

. JSCS 

30 

♦ I 

31 

♦ VOCG 

12 







. TEMPORARIES 











.EFF16 

70 












EFFIC ' i NO ERRORS detected ') 


CHAR LN'JGXB.FOR 


FORTRAN V.4(21B) /K1 


I 5-SKP-75 


9:D5 


PAGE 1 




SUBROUTINE CHARCE.T, ISC, VOC, FILL. EFF.ICOOE.VOCO, ISCO, FIFO! 



REAL IL, ILl , IMP.lHPl , ISC , M 1 , M2 . ISCO 

B0B03 


. FlFGtX, r)=(2.2784(X+4b00. )**-.112)-C6,666f>4)*T 

00004 


FV C X, y )=. 5476454 (X*4. 048 )-t 2. 45E-3) 4V 

00305 


FlCX.y 1=(( .14814Cy4».3631+.38814X 

00006 


FVS(X, y)=(2.9E-3)4( X4». 172)-f2.23E-3)4y 

0000-7 


FIS(E,n = (A0 + T4(Al +T4(A2 + T4(,A3+T4(A4+A54T) in 14E/10. 

00008 


___AQ=..8147 27, ; Al = . 1087 1 3E-2 : A2=-.695706E-5 

' 00009 


A l = .226603E-7 : A4=.17109E-9 : A5=-. 1 44039E-1 1 

■ 00010 


IFCICOOE .GT. 1) GO TO 22 

00011 


FTLL = FIF0.+,FIFGCE,T) 

00012 


voe=vbco+FV(E.T) 

00013 


ISC=ISC04F1(E.T) 

00014 


GO TO 21 

00015 

22 

VOC=VOCO+FVS(E,T) 

00016 


ISC=ISC0*F1S(E,T) 

00017 


FILL=(FIFOtFlFG(E.n 1 

00018 

23 

CONTINUE 

00019 


EFF=FILL*V0C4ISC4( ( 1 .0E61/E) 

00020 


RETURN 

00021 


ENli ■ 


SUBPROGRAMS CALLED 

FIFG 

FV FVS p-j--- -pjg 


SCALARS 

AND 

ARRAYS C "4" 

L HP. 

EXPLICIT DEFINITION 

- NOT 

RKFFRENCED 1 



4A4 

1 

%IL 


♦ FIFO 

2 

4EFF 3 

♦ T 

4 

.Fciwn 

■ .F0001 

6 

‘FILL 

7 

♦ A3 

10 

4E 11 

•F0010 

1 2 

UMP 

4IC0DE 

13 ■ 

■4 AO 

14 

%M2 


.F0007 15 


16 

.F0003 

4A2 

20 

ISCO 

21 

%M1 


ISC 22 

♦ A1 

2 3 

ifclMPl 

4VDC0 

24 ■ 

.F.0006 

25 

.F0004 

26 

•F0002 27 

•F0000 

30 

%ILl 

4V0C 

31 

4A5 

32 







TEMPORARIES 









.CHA16 

33 

.00001^ 

34 

•00001 

35 






CHAR C NO ERRORS DETECTED ') 



TNin V.4C21^) /Kl lS-SEPr75 9:05 PACK 1 


00301 

00307 

03003 

00004 

00005 
00008 
00007 
3000H 

00009 

00010 
0001 1 
00017 
03013 
00014 

00315 

00316 2 
00017. 

0001 P 

00019 72 

00073 
00021 


5UB90IJ n JC IfUTi (E, r, VOC. TSr.KF.eFP. ICODE. VOCO, ISCO.FIFO) 
9EAIj r^P, 1 SC» TL,M1 ,M7, ISCO, 11j1 . IMP I ,.M7N . «1 Pi 
KtF'^{ K. r ) = (7.?7«^( Kf 4 500, )**--l 1 2 1 - ( 6 . 666E-4 1 ♦! 

Ftf()^,O = ,&47 545*(X»*.04ftl-(2.45t-i)7y 

FLCX»t ) = U. U8 ) + (¥♦♦. 3b 3)f.3Hfl)+X 

FVS( X, Y )=( 2.PF- n 1 72 )-C 7. 2 3E-31 *Y ' 

FIS(E, r) = (AO+l*(AUrt(A2+r^(a34T*(AA + AS*T)))n*E/l0. 
AO=.9tt727 ; A 1 = , 1 ’013 7 1 3E-2 ; A?'=- , b9573hF-5 

a 3=, 2?bb'3 U:.-7 : A4=. ! 7n^9h>9 ■ ; A S = - . l 4 A039E- 1 1 

I SC = EFF4(E/ n .t'-Fbn /(FF# JOC) 

1K( rCOUF ,Gl’. n GO TO 2 
FlFO = F»-"-FlFG(E,T) 

VOCO=VOC-Fv(t,T) 

ISCO^ISC/*-'! (E.T) 

GO ro 22 

VOCO=i/OC-FV3CE,i) 

FlFO^Ff-FlFGU:,!) 

ISC0=£SC/FIS(E, D. ■ 

COJTI.JUr: 

3Et,NN 

EN') 


SUBPROGRAMS CAuLF.i} 

FJFG 

FV FVS FI FIS 


scalaps awo arrays I "*'* pjo explicit DEFiNnioN - hiOT referenced 1 


t A4 

1 

%lb 


*FIFO 

2 

♦ EFF 

3 

♦ T 

4 

.F031 1 

5 

-F3331 

6 

i;M2N 


♦ FF 

? 

♦ A 3 


♦ E 

1 1 

.F001V1 

12 

feIMP 


♦ ICOOE 

1 3 

♦ AO 

14 

1 N 


4;M2 


.F3007 

15 

. F3005 

t 6 

,F^003 

t 7 

♦ A2 

70 

[ SCO 

21 

IsMl 


ISC 

22 

*A1 

23 

4ilMPl 


♦ >/OCO 

24 

, F 

25 

«F00«)4 

26 

.F0002 

27 

.F00O0 

30 

[ LI 


♦ VOC 

31 

♦ A5 

32 






temporaries 

.lNJlb.33 .O00vlvj 34 .O0W01 35 


1 N 1 f I 


C 40 errors OEIECrED ] 



Oi 

Oi 


. MAXSI FORTRAN V.4f2lv») /Kt 15-SFP-75 9:9^b PACK 1 


3?>^33 

Ti^idTib 
7i^z?n ■ 


’sui^Kou r J -Mr: maxs r (f, ff'iax^temp, AR fCA . voro, tsco,fifo) 

r^FAb ISCO, 1 L, laC, 1 "P 

r=?s. ■; icooF =2 

CAbb CHAR(F,T,1SC, '/()C.Flbb,FFF, tCOnF. VOCO. T SCO, FIFO) 
FF^AX^FFF ; TFMP=?S. 

Rt TURN 
£'<0 


SUBPROGRAMS CAbbFO 
CHAR 


SCAbARS 

AMO 

ARRAYS [ MO 

KXPblCir nFKJMlTIOl'i 

- "%** NOT 

RtFFRFNCFD I 



lAREA 


fb 

♦FIFO 1 

♦ KFF 2 

♦ r . 

3 

♦ FILL 

*F. 

5 

♦(COOK h 

^ l -iP 

♦FFMAX 7 

ISCO 

10 

tsc 

*voco 

12 

13 

♦VOC 14 






TKMPORAR [P:S 
. .M4X1 h 1*3 

MAXSI I ERRORS OtirECfFO ) 



ONTWO [iN3GXB.F0R 


FORTRAN V. 4(210) /KI 


15-SEP-75 


9:05 


PAGE 1 


Oi 

-3 


00301 . 
30302 
00303 C 
30304 C 

00305 

00306 

00307 
0030B 

30309 

30310 


ggaii 


.30312 

03313 

30014 


30315 

33016 

00317 


30318 

00319 

30320 


30321 

30322 
00323 


30324 

00325 

00326 


00027 

30328 

00329 


00330 

30331 

30332 


30333 

00334 

00035 


30036 

00037 

00338 


00039 

00340 

00341 


00342 

30343 

00344 

30345 

00346 

00347 


.30348 

30349 

00350 

00351 

00352 
00?53 


00354 


101 


_2 

102 

_3 

105 


. 77 


4 


.5 



SUBROUTINE 0NT40(.TMAX,CFPSI , IK , I N , FPS . VOCO. ISCO. FIFO. HI , TF, ETA3M, 
2rEMP .AREA) 

IF IK = 1 THISj silRROUriNE WORKS ON SYSTEM 1 
IF IK = 2 this' SUBROUTINE WORKS ON SYSTEM 2 
REAL IN. ISCO J HP , I SC , I L 

FCI) = ’(1 .-E>F)4'en-H 1*(T-TAIRK)-EPS4.5IGMA*( n *»4)-(TAIRK»»4) ) 
DELT=(TMAX-TFD/10. ; S IG*! A = 5 . 6697E-8 : JAIRC=25. 

_OEL_T.= 5_. 

EN=IN»,84 

TAIRK=TAIRC+273. ; TFK=TF+273. 

IFCIK .G f. n M TO 2 
WRITE(3,101) , 

FORMAT!//, T30, 'SYSIEH I ' , / , T 10 , ' FLAT PLATE COLLECTOR FOR THERMAL F. 
.2NER.G!f, A.NO' ,/,T10, 'FLAT PLATE COLLECTOR (SILICONJ FOR ELECTRICAL EN 
3ERGY’ ,//) 

GO TO 3 

WRITE(3.1 021 

FORMAr(//,T30, 'SYSTEM 1 1 ' , / , 110 . ' FLAT PLATE COLLECTOR FOR THERMAL 
2ENERGY ',/, ri0, 'HAVING SILICON (EFFIC. FCN. OF TEMPERATURE)',//) 
_j:_P'NTINU£ 

WRifE('3, 105) 

FORMAT (T3, ' TEMP 'T9, ' R ATE/ AL ' T20 , 'QTHERM'T31 , 'QT0TAL'T42, ' FRAC'T48, 
2'FRAC'I57 , 'F'T64, '«U/AL'T7b, ' 0 TOT ' T86 . ' FRAC T93 , ' FRAC ' T101 , ' FI ' ,/, 
3r3,'(C)'T9,'W/M2.K'T20,'W/M2'T31, ' */M2 'T57 , • s/w • T64, ' W/H2 ' T7h , ' w/M 
42 ' T86, 'ELECT42, 'ELEC T4 8 , ' THERM 'T9 3, ' THERM • T100, ' S/W ' T1 07 , ' EFF' , 
51111, 'F3($/H)' ,//) 

00, 60“l=i , 10 
TK=TF+IAD£LT+273. 

TC=TK-273 . ' 

1F(IK ,GT. 1) GO TO 4 
IF(ETA3M .NE. 0.) GO TO 77 

CALL M AXSI(1N,EFMAX,TEMP,AREA,V0C0, ISCO, FIFO) 

ETA3M=EFMAX ; EFF=EFMAX 
BBR=( rK*»2)+(TAIRK4»2 ) 

AAA=BBB *((TK*42)-(TAIRK4*2) ) 

JUAL=EN-H14(TK-TAiRK)-EPS4SIGMA»AAA 

3ELEC=ETA3M»EN/100. 

_go1to_5 

CALL CHAR(EN,TC,ISC, VOC,FF,EFF,2,VOC0.1SCO,FIFO) 

EFF=EFF/100. 

Q UAL=F ( TK ) 

3ELEC=EFF4EN 

EFF=EFF»100. 

_IF(_QUAL .LT. 0.) GO TO 50 
3T0T=QELEC+'0UAL 

PE=OELEC/QrOT : PT=OUAL/aiOT 

RAT E=QI)AL/(TK-TFK) 

ETA2=aUAL'/E,N 
CFP=400.4ETA2 +10. 

3rHEHM=.54(l.-TFK/rK)W0UAL 
' 3T6TAL=0'E'LEC+QrHERM 
PQE=3ELEC/0T0TAL ; PTH=3THERM/0T0TAL 
Jl^CCFP+^FPSDyOIOTAL 
F2=(CFP+CFPSI)/QT0T 
PE=PE4100. ; PT=PT4100. 

POE=PQE*100, ; PTH=PTH»100. 


00355 

30356 


05 

00 


O'lirwO LN3GXB.F0R FORTRAN 4(213) /KI 15-SEP-75 9:05 PAGE 1-1 


00357 F3=(CFP+CFPSI)/0ELEC 

0005R .»(R.IJE.C3_,J04JTC,RATE,OrHERN, 0T0TA[,.P()E.PTh,F1 .DUAL, OTOT , PE , PT , 


00059 


2F2,EFF ,F3 


00060 

104 

F0RMAT(1X,F5.1 .1X,F10.4,1X, 1.PE10 

. 3.tX.Et0.3,lX.0PF5.2.1X,F5.?,lX 

00361 


2F9.3.1X.1PE10.3. 1X.E10.3. 1X.0PF5 

.2, 1X.F5.2. 1X.1PE10.3. IX,0PF5.2, 

00062 


31X,F10.2) 


00363 

50 

CONTINUE 


00364 


IFIIK ,E^ 1) «HITE(3,80) EFMAX, 

TKi'lP 

30365 

80 

FORMATCllX, 'MAXIMUM EFFICIEMC7= 

'.Ft0.3.* AT TFMP. OF *.F10.3* 

00366 


2' DEGREES CENTIGRADE') 


• 00367 


XRITEUj 107J OELEC.CFPSI. 


30368 

107 

FORMATdX, 'QELEC = ' , 1 PE 1 0 . 3 , 2X , 

' «/M2 COST OF SILICON = 

00069 


2 0PF10.2, ' S/M2' ,//) 


00070 


aSTURN 


00371 


END 



subprograms (TAUUEO 


CHAR 

MAXSl F 

SCALARS .AND ARRAV S [ S0_ EXPLICIT DEFINITION - '•»» NOT REFERENCED ) 


♦CFPSI 

1 

♦ AREA 

2 

%IL 


♦ FIFO 

3 

♦ EN 

4 

♦ EFF 

5 

♦ POE 

6 

iFOTOT 

7 


10 

»FF 

1 1 

♦ Ml 

12 

♦ OUAL 

1 i 

♦ IS 

14 

♦ PT 

15 

%[MP 


♦ TFK 

\h 

♦ DELT 

17 

»TF 

20 

♦EFHAX 

21 

♦ PTH 

22 

♦QTOTAL 

23 

♦ TAIRC 

24 

ISCO 

25 

♦ ErA2 

76 

*TK 

27 

♦QTHERM 


♦ F3 

31 

♦ CFP 

32 

♦OELEC 

33 

ISC 

34 

.50000 

35 

♦ F2 

36 

♦ RB6 

37 

♦ voco 

40 

IN 

41 

♦FTA3M 

42 

♦ TEMP 

43 

■♦RATE ’ 

44 

♦ TMAX 

45 

♦ voc 

46 

♦ I 

47 

♦ FI 

5M 

♦ TC . 

51 

♦ PE 

52 

♦TAIRK 

53 

♦ EPS 

54 

♦ AAA 

55 

♦ SICtMA 

5f> 


TEMPORARIES 

.ONTl'6 325 

ONTWO r NO ERRdRS O'^TECTED ] 



THRFO liN.JoXrt.FOR FORTRAN V. 4(210) /Kl l5-8FP-^‘75 9:05 PAGE 1 


Oi 

CD 


7i?idd\ 

TiVtZdi 

ZTiZZb 

ZZZZh 

ZdZZl 

ZZZZB 

z?zd9 

ZdZ] 1 
Zd'd\? 
3 

ddd\^ 

5 

?ZZ) 6 
ZZd.\ 7 
dZAXH 
ZdZX 9 
ZZZ7Z 
7)0321 
33322 
3032'3 

33324 

33325 
333?f> 
33327 
3332B 
33329 
33333 
33331 
333.32 

33333 

33334 

33335 
33335 

33337 

33338 

33339 
33343 

33341 

33342 

33343 

33344 

33345 

33346 

33347 
3334H 
33349 
33353 

33351 

33352 

33353 

33354 

33355 

33356 


SUBROu riNE THRFOC IK , TMAX pCFPSl , I N , FPS , VOCO . ISCO . FI FO . H 1 ,TF,ETA1 , 

2 A R E A ) 

C [F 1K=^ THIS SUHROUriNE -^ORKS ON SYSTEM 3 

C .tE IK=1 THIS SUBKO'JVTNE ^lORKS OW SYSTtM 4 

, real In ; 1 SCO. I MP. IL, ISC 

S( D=( 1 . -EFE )»E-2. »H1 ♦( I’-TA 1 PK > *2 , *F P.S ♦S JHMA * ( (T*»4)-(TATRK^»4))- 
Pl = 3. 1 41 : SIGMA = 5.b69 7F:-B : TATrtC=?5. 

OELV=( r-fAX-FEJ/^O. : TA lRK = TA-[PC+2 73, : rFK = TF + 273.. 

OEL(’ = M, 

TFUK ,GT. ? ) GO 10 1 
'VRI TE( 1,133) 

133 FOH-vArc//. i'3o, *srsrF:^ in ' ./. [lo. tumcentrator with ga-as solar ce 

2LLS' ,//) 

JO ro 2 

1 ^R I TE( 3 . 131 ) 

131 FORMATC//, no. 'SYSrE'^l 1 V* . / » T 1 3 . ' CONC EN I'R ATOR WITH SILICON SOLAR C 
2ELLS'.//) 

2 COntNUE 

)0i=i3. : crH=-23. 

00 52 J = ni5 
CrH=CTri+001 

1E(J .GT. 4) Ci'H=crH + 20. 

1E(TJ .or, H) .AND. (IK ,E0. 37) CTHs 30P- ( J-8 ) * 1 00. 

lECtJ.Gl'.RI .AND, (IK ,E0. 4)) Gu TO 52 

CG=CEPSl 

CC0N = M3. *F,lAi *R0. 

wR rVE( 1. 132. ) CTri, ETAl.CCntM.CG 

132 F'JR'^AT(/.U. 'CONCti'jrRATlON - E 1 13 . I . 5X ,* FURNACE FACTOR = ETAl = 

2 ' ,E1 3.4, /. 1 X . 'COST OF CO^'CEfJ TRATOR = *, F 1 0 . 2 , 2X , * COST OF SOLAR CEL 
2L = 3 ’.l-’n.-;,//) 
iRI IE( 1,134) 

134 FOR -1 A T( n, 'TEMP' 19, ' RATE/ AL' 120. ' UFLEC'TIl , *QTHERM'T41 . 'OTOTAL 'T5 3 
2 . ' FRaC T59 F RAC T6b EFF' 172 F'TP 3. '0U/AL*T94, 'OTOT'T104, ' FI 'T1 1 
34, ' FRAC ri2l . ' FRAC , / .T 1 , ' 1 C) ' T9, ' vJ/ft2.K 'T20. * VI/M2 ' T 31 , ' W/M 2 ' T4 1‘, 

4 ' W/M2 ' f 5 i ELEC f 59 , ' THERM ' T7 2 , ' S / w ' Tfi 3 . * W/M2 * T94 , ' N/M2 * T1 34 , 
5'5/K'^'rm,' FLEC fr2l.'THERM'Tl27,'FUS/wn,//) 

E = ErAl 

DO 52 1=1,23 

rK = TF+l*DEL J +271. 

rC=TK-2 1 1. 

1F(IK .GT. 3) GO TO 5 

CALL CHAR( E:,[C, ISC, \/0C. FF, EEF, 1 , ^/OCO. I SCO. FIFO) 

JO ['0 6 

5 CAIjL CHAR(E,TC, ISC, V0C,FF,EFF.2, VOCO.-ISCO.FIFO) 

6 OELEC = EFF + E/( ly'3.*CrH) 

■:ef=eff/-io3. 

)UAL=G(1’K) 

IF()UAL .Lr. 0 . ) GO ro 52 
RArE=OUAL/( I’K-rFK) 

C2 = C«:i/Cl'H 

0THEK3=,5E( 1 .'TFK/ TK ) +QUAL/CTH 
, J2=0JAL/CTH 
0T0T = :r2+ JELKC 
jrOPAL=QElEC+OTHERM 
■PE=0ELEC/0 TOT ; 91 = 02/0701 

?OS = OF:L£C/OTOTAL ; PTHsOIHERM/OTOTAL 



THRFO 


f'OPTRAN V.4C>1{^) /Kl 


1 5-SRP-75 


9:05 


PACK 1-1 


FOR 


!1035H 

dddhl/i 

2)0?161 

dd3b3 

0?ia6fl 

00?)69 

00070 


F=(CC0'>JfC2 )/oroiAi.. 

F2=CCCO*^+C2 )/OTOT 
F2=F2*1000, 

PR = PF^MO. : PT = pr*100, ; EFK = f:FP* 1 00. 

P l|L = PO^ + t00, ; PTrisPTHflvlw. 

r 1=lCCOMtC2)/Oei.EC 

I I'EC 1 » I )TC, RATE. OELEC. OTHER M,0|'Or A I.. POE,PTH.KFP.K.QUAI..QTOT*K2 
2,PE,P[' .FJ 

10i FOR‘-iA rc 1 < . P’5, t , IX. 1 PE10. i . 1 X .Elrt. 3. 1 X . El 0. 3 , t X . P' 1 0 . 3 . 1 X . 0PF5 . 2 . 1 X . 
2F5 . 2 , 1 X , F5 . 2 . I X . 1 PF 1 0 , 3 . 1 X . P: 1 0 . 3 . 1 X . F 1 0. 3 . t X .0PF1 0. 2 .IX. 

3 1 X , K5, ^ . 1 X ,F5. 2, IX . F7 .2 ) 

52 CONTI OiJiC 
REI'URJ 
EMO 


SUBPROGRAMS CAUlPM 
CHAR 


SCALARS 

AJ ) 

ARRAYS 1 

NO 

?XPL1CU' OKKIUITU'W 

- 

NOT 1 

REFERENCEn 1 




♦CFPSl 

I 

%AREA 


■*IL 


♦ FIFO 

2 ' 

♦ EFF 

.3 

♦ PQE 

4 

♦ QTOr 

S 

, K'0000 

6 

♦ FF 

7 

♦ HI 

10 

♦ PI 

11 

♦ E 

12 

♦ OUAf, 

1 3 

f IK 

1 4 

*p r 

1 5 

%IMP 


♦ TKK 

16 

♦ OELT 

17 

*TF 

20 

♦ P PH 

2 { 

♦ 0 rorab 

22 

♦TA IRC 

23 

ISCO 

74 

*TK 

25 

♦ OTHER'* 

2b 

♦ F3 

27 

*j 

30 

♦OELFC 

31 

JSC 

32 

.S«O0l 

33 

♦ ETAl 

34 

, S 0 0 -0 VJ 

35 

*F2 

3h 

• voco 

37 

IN 

40 

♦ C2 

41 

♦ RATE 

4 2 

♦CCON 

4.3 

♦ TMAX 

44 

♦ VOC 

45 

♦ CTH 

46 

♦ I 

47 

*F 

50 

+ rc 

5 1 

♦ PE 

52 

♦ lAIRK 

53 

♦ ODl 

54 

♦ EPS 

55 

*02 

56 

*SIG«1A 

S7 

♦ CG 

bO 








lE'^POPARlES 


.I’HRIB 302 


THRFO 


I \)j ERRORS OFTECrFO ) 



FIVE LN3GXB.F0R FORTRAN V.4C2ia) /KI 15-SEP-75 R:0b PAGE 1 


80201 SUBROUTINE FIVE! TM AX , IF , I X , CFPSl , IN , EPS, VOCO. XSCO . FIE'O, H 1 .ETAIN, 

00002 lyAl., TEWP.AREA) 

00003 C THIS SUBROUTINE WORKS ON SYSTEM S 

00004 real, in; isco, ISC, imp, IL 

00005 GlTJ=E-2,»H14(T-TAIRK)-?.YEPS*SIGMA»f (T»»4)-tTAIPK* + 4)) 

00006 NRITEIJ, 101) 

00007 , 301 FORMAT! //, T30 ,' SYSTEM V ',// ,Tt0, 'CONCENTRATOR FOR THERMAE ENERGY 

0000 8 .2., AND SEPARATE SILICON SYSTEM FUR ELEC. ENERGY') 

00009 PI = 3. 1415926536 : S IGMA = 5 . 6697E-8 : TAr.RC=25. 

00010 TAlRK = TAlRC4'27 3, 

00011 ^rFKj=:rF.t_27 3 , . 

00012 DELT=CTMAX-TF)/1 5. 

00013 001=30. ; CTH=-20. 


..OQ 5.0 0=1,8 


00015 CTH=CTH+u01 

00016 IF(J .GT. 4) CTH=CTH+20. 

00017 CC()N = 1 0,0^. + ETA 1 ♦ 8 0 . 

00018 CFP=CFPSI 

00019 IFIETA3.M .GT. 0.) GO TO 3 

00^20 CALL MAXSHIN.EFMAX, TEMP, AREA, VOCO.ISCO. FIFO) 

00021 ErA3M=EFMAX 

00022 3 OELRC=ETA3M+(.S4)»)N/100. 

00023 WRiTEC 3,102) CTH,ETA1 , CCON , CFP , ETA3M .UFLEC 

00024 102 FORMAT!/, 1 X, 'CONCENTRATION = ', FI P. 2 , 2X ,' FURNACE: FACTOR : ', 

00025 2F10.3,/,IX, 'COST OF CONCEN rRATOR = ' . F10 . 2 , 2X , 'COST OF SILICON = ' 

0M26 3,F10.2,2X, 'MAX EFF OF SILICON = ' , E 10. 2 , / , 1 X , 'OEIEC = ',fI0.2.//) 

00027 ' WRITE! 3, 104) 

00028 104 FORMAT! T3, 'TEMP 'T9, ' RATE/AL'T20, '02 ' T31 , ' OTHERM ■ I'41 , ' QTOTAl, ' 1 53 

0 0029 ' _2,.1FRAC'T59, 'FRAC"rb6, ' EFF ' T72 , ' F ■ T83 . • OU/AI,' T94, • OTOT ' T104 , ’ FI 'Til 

00030 34, 'FRAC'''t 1'21 , 'FRAC ,/, T3, ' ! C ) ' T9 , ' W/M2 .6 ' T20 , ' w/M2 ' T31 , 'w/M2'T41 , 

00031 4'W/M2' 153, 'ELEC T59, ' THERM '. 1 7 2 , ' S/w ' T8 3 , 'k/M2'T94, 'w/H2'T104, 

- _ 5J S/Kw' n 14, 'ELEC'Tl 21 , ' THERM ' T 1 27 ,' F3 ! S/Wl '.// ) 


00033 


00034 

DO 1=1,15 

00035 

rK=TF+I40ELl+273. 

0003.6 

TC=TK-27.3. 

00037 

3UAL=G!TK) 

00038 

IF!OUAL .LV. 0.) GO I'O 50 

00039 

RATE=QUAL/!TK=TFK) 

00040 

3THERM=.5«! 1 .-TFK/TK )tOUAL/CTH 

00041 

OTOTAL=QELEC+OTHERM 

00042 

PaE=OELEC/OTOTAL ! PTH=()THERM/OTOTAL 

00043 

F=(CCON+CFP)/OTOTAb 

00044 

32=QUAL/CTH 

00045 

3r0T=02 + (iELEC 

00046 

PE=OELEC/OTOT ; PT=Q2/QT0r 

00047 

F2=!CC0N+CFP)/0T0T 


00048 PE=P£*100, ; PT=Pr*100. ! EFF=EfE»t00. 

00049 POE=PQE4100. : PTH=PTH»100. 

000_50 F2=F2410'00. 

00051 EFF=ETA3M 

00052 F3=!CC0N+CFP)/0ELEC 

00053 WR1TEI3,103)TC,RATE,Q2 , OTHERM ,0T0TAL. POE, PTII , EFF, F,OU AL , 0 rOT , E2 

00054 ^ ' 2;PE,PT ,F3 

00055 103 FORMATI1X,F5.1,1X,1PE10,3,1X,E1H.3,1X.E10.3,1X,E1».3,1X,0PF5.2,1X, 

00056_ ■ 2F5.2,IX,F5.2ilX,lPE10.3,lX,E10.3,lX,E10.3. 1X,0PF8.2,2X, 



FIVE LNaGXS.FOR FORTRAN V. 4(21511 /Kt 15-SER-75 OtOb PAGF 1-1 


00357 31X,F5.2.1X,F5.2,1X,F7.2) 

00058 50_ CONTINUE 

00059 RETURN 

00060 end 


SUBPROGRAMS CALiLEO 

_ ^ _ 

MAXSI 

SCALARS AND ARRAVS [ NO EXPLICIT DEFINITION - •'%» NO! REFERENCED' 1 



*arr:a 

2 

%IL 


♦ FIFO 

3 

♦RFI* 

4 

♦ POE 

5 

♦OIOT 6 


7 

♦ HI 


♦ PI 

11 

♦ E 

12 

♦ OIJAL 

1 3 

tIK 

♦ pr 

14 



♦ TFK 

15 

♦ DELT 

1 6 

♦ TF 

1 7 

♦ EFWAX 7.0 .. 

♦ PTH 

21 

♦OTOTAL 

22 

♦ TATRC 

23 

ISCO 

2.4 

♦ TK 

2S 

VQTHERM 26 

♦ F3 

27 

♦J ■ 


♦ CFP 

31 

♦OF.i.t:c 

.32 

%ISC 


.S0001 33 

♦ ecAi 

34 


35 

♦ F2 

36 

♦ voco 

37 

TN 

4W 

tETA3M 41 

♦ TEMP ■ 

42 

♦ RA TE 

43 

♦ CrON 

44 

♦ TMAX 

4S 

♦ CTH 

48 

»I 47 

♦ F 

50 

♦ rc • 

SI 

♦ PR 

52 

♦TAIRK 

53 

♦ f>I)l 

54 

♦EPS 55 

*Q2 

■56 

♦SIGMA 

S7 








TEXPORARI.ES 

_t.LU.16_3.® 2. 


FIVE 


( NO ERRORS OETECTEU 1 



-Cl 

GO 


MAIN. LNjn^-^.KOW KOKT«A'J V.4(9l^) /K| l^-Se■P-75 9:i’S PAGKl 


7t'A?i?'\ C Sl^SrP:^ AMALY«iIS 

TiTiTid'l ■)Tn:<S10'l CHPk( 2M , VO[,l 2b 1 

?i?>7i'A-i kf-;Alj ISC, ISCO, u,. [SCI 

ktAij ISCOG, iSCOS,I^J 
kt-Al/ l.SC".,lbCS 

7t?)'A'Ab J\ = 2o ; ; h,rAiM = t^. ; lTl^*F-.=t^. 

?A??>7 akK’.A^S ■)'/>. 

?5p?^H c ^‘•:ao S()lak CKbu oarfl 

7/^/1 ti ['-'i*:=m‘iF:M 

A?>AUl -ItAlU 2, I , VOC, ISC. F K, tf'K. AM(w imDE.LAST 

7 AA r'Ok -u rc / (F , ?( rs ) ) 

2 1 l>J( i. I . VJC. ISC.FP.KKF, A-0. 1 rOPF.I.AST 

7i‘A?i] ^ 2-?1 POk'M rt I .< . 7( FI 2,^) . 2 t I S n 

CAijb IM I I’UK. r, Vf.)C. ISC. FK, CPF. ICOOF. vocn, IPCO.FIPO) 

?3^JS rpIlCuPF ,or, 1> CO fn 7^1 

'A'A?>]h rcc 3. 202 1 

7iA?\'l 202 !'l // , U . 'GA-AS SOLAP CFCL DATA'.//) 

Tiooif* i/;jCOG=vnC(i ; rsCOG=ISCO ; FIKuG=FlKn rARfAG=AREA 

{ )’F( 3 , 2H1K. r. VOC. ISC. FF. CPF, VOCO. ISCO.KIFO 

,i?a020 2o3 FOkoA rt 1 X, ' I I-I VFijStrV - ' .FM.b. IX, 'TFVP = * ,F'1H,2,2X. 'VOC =',KU1.5 

Z7>A?] 2,2<.'1SC = ' .K1 0,b, / . 1 X , ' KI LI, FACTOH = * . F 10 . S , 2X , ' EKK TC f P-^CY = 

?>??22 3F1/J.S.2X, / . IX . ' voro, TSCO.KlPn = ' ,.3f 2XF1 4,f») ) 

3 00 2 3 G 0 F 1 > 7 o 2 

33024 7 01 I r.K( 3 , 204 ) 

33?2b 204 FOkU n //, 1. X , ' SIM CON SOLAK CFliL HATA'.//) 

333?f> VOCOS = VOCO ; rSCOS=lSCO ; FiFOS=PrFLt fARt:AS=APFA 

333'2 7 I. t'FC 3 . 2 ).3 )K, r. VOC. LSC. PF . PFF. VOCO . 1 SCO , K I FO 

3-332H 702 lF(iri^lF ,L1. 2) GO TO 7‘0<0 

33329 CAuL FKFLC( VOCOG. I SCOG , F I FOG . VOCOS , 1 SC(^S . F IFOS ) 

33333 C 4FA0 SYSTF'*'- PARAMprCHS 

33331 C rp= I'FopF^a 1 iiwj/ op FLUID 

33332 C FI A1=F‘JW''jACk f aCTOk - ,i AMD .7S 

33333 C V = COST OF (;A-AS (F + 34 DOLLARS) V = 90 . »0 . 70 . 60 . 50 . 0« 40 

33330 C I V = l -Jl’KM.Sl 1 V OF SOLAR KAorariOU lo w/m2 

3333S C FPS = F OM I SsS 1 r [ V I n OP comcfni’rator 

3333b C 91= CO gVFCT'lOi'J HEAT LOSS COfJFFJCIFMT In H/M2.K 

33337 99 REAO(2,9JM TF , E TA 1 , V . EPS , H 1 . IN .LAST 

33338 931 FOR A T f h C P 1 0 . 0 ) , I 5 ) 

333 39 ^R FTFC 3 , 9H I )TF . E FA 1 . V , EPS , H 1 . 1 v , LAST 

33340 9H1 FOR'^A rC 1 < .b(E1 2.5) . I S)' 

3334] IFCrF .Gr. Hu.) GO TO 932 

33342 JJ=l 

33343 GO lO 434 

33340 9 32 J.l=3 

33345 934 IK=JJ 

3334b C COS = Cost of SILJCOM 
33347 COS=10O0, 

3334H CSPFP^CDS+54. 

33349 IPCIK ,GT. 2) GO TO 9)6 

30353 call 04 PRO ( t OO . . CSPFP , 1 , I N , tPS . VOCOS . I SCOS . FI FOS . H 1 . TF . ETA )M. 

33361 2 rE'^P, AREAS) 

33352 call 0'U>jO( 1 03. .COS , 2 , i N . P P-S , VOCOS . 1 SCOS , F IFOS , H I ,TF,ErA3M, 

33353 2 TE'Xp. AREAS) 

33354 CALL THRFOI 3, )00. . ( V* 1 ,0E4) , L N , KPS . VOCOG , 1 SCOG . F T FOG • H 1 , rF.FTAl , 

33355 2AREAG) 

33356 CALL I'HRFO ( 4 . 1 93 . . COsS , I M . EPS . VuCOS , I SCOS . F IFOS, H t .TF .ETA1. 


-Cl 


MAIN, L^4)GXR,F0W fORTKA'g /KI 1 5-SFP-75 9t«5 PAGE 1-1 


03358 

33359 

33363 

33361 

33362 

33363 

33364 
33366 

33366 
33^*67 

33368 

33369 
33371? 
3337 1 

33372 

33373 

33374 

33375 

33376 

33377 

33378 

33370 
33383 

33381 

33382 
333H3 
333H4 
333R5 
33086 

33367 

3 3 3 P 6 

3 3 3 M 9 
3339? 


2AR-:a.5) 

CAIjlj F WK( 6Hvl, ,TK. 6,CSPF‘P , I N , KPv5 . VOCOS . I SCOS . FIFOS . H 1 ,ETA3M.FTA1 , 
2 r F "1 P , A F F A S ) 

GJ i'O ^52 
9^6 CO>lTi:H)F 

call r-lRFOt 3 , .^/O, . ( V* I ,3K4) . J m , RPS . V OCOG • !SCOG.Fi6'OG,Hl . rF.ETAl , 

2ARii;At;i 

CALL r H F FO c 4 , 1 ri (3 . . COS , f , F PS , V OCUS , I SCOS , F T FOS . H 1 • TF . E T A 1 , 

2 A P K A S ) 

CALL F t \/E( 5^1?. , i’F, S.CSPFP. i »•» .EPS . VOCOS. I SCOS , F I FOS . H 1 .6:TA3M,FTA1 , 
2l’F. -IP.AHFA.S) 

bS2 IFCLASr ,Gl’. /) GO TO 99 

LF(L.ASr .GP. -4) GO J’O 631. 

0.) 6.1? KK= 1 . 7 
FFF=W. + (<><-l )*2. 

F=10vM, J T=26. ? V0C=1..O : 1SC = V>.W19 J FFsO.flS : AMO=JA. 

rcOOL’ = 1 

CALL i 'I r n (K/r. f sc,ff,fkf\ icoof , voco.isco.ftfo) 

.>p t IK( 3 , 2 12) 

;:)COGsV0CO ; ISCOG=USC0 ; FlFi)G=FlFO JARFAG = AREA 

^■RI lr( 3,2/;3 IF. r,vOC. I.SC.FF.FFF. VOCO . 1 SCO, FI FO 

CALL rRPr 0( .3 . 33^. , ( V ♦ 1 .^F4 I , ! .V . F.P>S , VOCOG . TSCOG. Ft FOG . HI . TF.FTAl , 



2 A rt F. A G ) 



532 

COMTl J'JF 

CALL I 'i I 1’ 1 ( 1'3'0 0. . 25. . 1 .3, LSC. .H6 

.16,.] 

1 ,VOCOG,;srOG,FIFOG) 

531 

CO i 1 1 •'J )K 



63 3 

RKAU( 2 . 9 31 ) PF. MJI'A I . V.FPS. H 1 . T.iv, 

liAST 



C 0 A = 2 3 0 . 

CAt.L PHKF'U 3 . 333. , c V* I ) . I N . FPS . VOCOG , ISCOG , F IFOG , H 1 , I F.ETAl . 

2AKFAG) 

CAfOi IMRt-Ot 4 , M3. .CC'S, I •'J .EPS, VOCnS. I SCOS . F I FOS . HI .TF.FTAl . ARFAS) 
IFCLASr ,F:j,-1 ) GO VO 633 
STOP ‘ "'OF-iAlj FVO' 




F 1 y ' 

FFFl C PH 2F J I •'i 1 I’ I ON I '*JO 


SCALARS 

A '10 

ARRAtS t 

' -'in 

EXPLlCtP OKFINITIOW 

_ " S- '* 1 

k=OT 

reffremcfo 1 




♦vDCns 

1 

♦ AREA 

; 

ISCOG 

3 

♦ FI FO 

4 

ttL 


♦ JJ 

5 

♦ EFF 

6 

♦ r 

7 

♦ ff 

1 3 

♦ HI 

1 1 

♦ITIME 

12 

♦ f: 

1 3 

*TK 

i 4 

* 1 C 0 1.1 r. 

1 6 

%. i -IP 


♦VOCOG 

16 

♦ LAST 

17 

♦ AREAS 

20 

?;VOL 


♦ CSP 

21 

♦ PF 

2 2 

♦FIFOS 

23 

♦ V 

24 

ISCO 

25 

♦ K< 

2 6 

ISC 

?7 

♦ EPAl 

33 

♦apfag 

31 

.SvliAOW 

12 

%ISCG 


*F). FOG 

3 3 

♦ voc) 

34 

♦ AMO 

3S 

Ifj 

36 

I sens 

37 

♦ ETA3M 

40 

♦ COS 

41 

♦ VK-'P 

42 

4;1SCS 


♦ voc 

43 

%CURK 


♦ KPS 

44 

?:ISC1 


♦ 'V 9 

46 











APPENDIX B 


COMPUTER PROGRAM FOR ANNUAL ENERGY PRODUCTION 
OF SYSTEMS III, IV 



-3 

o> 


HMH," lNaiB2,rOP ' rORIRAN V.4C210') /K1 18-AUC-T8 


19122 PACE 


•eeei subroutine blocki 

eeeei dimension xk(S) 

_ 00001. COMMON ErriC0N,CELLl,.CElil2.H»Xl,XMiEIAl,XK,Tr 

00E04 NR1TE(3«233> 

00005 233 rORMAT{*t'ilX«'6A>AS SOLAR CELLS'.//) 

_000O6 NPITE(3i200) ErriC0N.CELLl.CELL2«EIAliir_ 

0000) 200 rORHAT(lX,'EFfICIENCT ■ * ,fS , 3, 2X, 'CONC COST ■'.Pl0.2>/i 

00008 21X.'C0ST CA>AS «' « ri2 . 2> 2X> 'COST SILICON «'ri2.2i/» 

_ 00009 JlX.'ETAl • 'ri0,2i2X, 'FLUID TEMP a'.riO.T.//) 

00010 CALL PCRF0R(l.ErF,C0N,CELLl,CELL2.N,XI«XM,CTAt>XK,TF) 

00011 RETURN 

_00012 END_ 


..COMMON blocks - . 

/,C0MM./(*2t) 

_irr _*0 ^CON 01 CELL! >2 CELL2. 01 N ♦« 

XI *S XH «« ETAl 07 XK 010 TP 020 

_.8UBPF06RAM.a_CAUED 

PERFCR 

■“SCAlTRa'AND ARRAYS t NO^XPLIClt OEflNirioiT • "%• No'rBWER'fNCED”) 


TEMPCRARIES 

“TBioTi' 46 

ilOCK.l . t NO ERRORS J)ETECIEP_1 



MAIN, lNSlB2,rOR 


rORTRAN V, 4(210) /KI 


18>AUC>TB 


19123 PACE 1 


00001 SUBROUTINE BL0CK2 

00002 DIMENSION XK(S) 

00001 COMMON Err,CON,CElLl,CELL2»N,XI,XH,ETAlfXK,ir 

00004 KRITE0.234) 

00005 234 rORHATC'l'.lXf 'SILICON SOLAR CELLS',//) 

00004 KRITE(3,200) EFr,CON,CELLl,CELl2,EIAl,TE 

00007 200 rORMAItlX.'ErriCIENCX « ',E8,3,2X,'CONC COST •',ri0,3,/> 

00008 21X,'C06T GA«AS «' ,F12 .2 , 2X, 'COST SILICON ■'ri2,2«/» 

00009 IlX.'ETAl « 'ri0,2,2X,'FLUIO TEMP ■',F10.2,//) 

00010 CALL PERrOR(2,Err,CON,CELLl,CELL2,N,Xl,XM,ETAt,XK(in 

00011 RETURN 

00012 END 


. COMMON BLOCRS . . _ _ . . . 

/,C0MH,/(+21) 

trr_ *0 CON ♦! CELLl__+2 CELL2 ♦! _ » _ _ +4 

XI 45 XM 46 ETAl 47 XK 4)0 TF 430 

subprograms called 

■PERFCR 

SCALARl and TRRAYsX"** Nb~EXPLiCIT omKITION NOT RErEREHCED' J 


TEHPCRARIES 
"TBLOIS 44 

BLOCK? . t NO ERRORS DETECTED) 



-J 

00 


WSrCF ' lN5LB?,rOR" rORIPAM V,4tai0) /KI 18«AUO»7S 19122 PACEl 


0BBB1 SUBROUIIHE mrOR’tlCbDE.trr ,CON,CEtil ,CE1-19 ,«iXi»X«.E2Al,XK,Tr ) 

00002 DIMENSION XK(I ) ,OAIM 1 1 2) • XIN( 1 2) , TA1R(12) f QE( 12) >QUSE( 12} rawORKC 12 

00003 2) 

00004 r(TC»X,Y»Z)«EIAR*(l..BETA*CTC*TRK))*ETAl#X*CTH-ETAl«X*CTH+ 2.« 

00005 2 HU(TC«Y)*2.*CPS»S1CHA*((TC««4)*CY««4))«Z*CTC-TFK) 

_00006 rP(TC,X,Y,2.)f?ETAR#BETA»ETAl«X»CIHf«,»EPa»31CHA*(TC**31 

00007 2 *2,«H14Z 

00001 DAYMCn»Jl, I DAYH(2>«20. ( DAYM())b 31. rDAYM(4)»J0, 

.00009 DAYM(5)»31, | DAYHt6)»30. I DAYMC7)«31, I DAYMt8)«31, _ . 

000iB DAYM(9]«30, I DAYMCl0)«3l, t DAYHai)a30. i DAYHC12)a31. 

00011 IAIR(1)>10, 1 TAIR(2)«l0i I TAIR(3 )b20i I TAIR(4)>20, 

_00012 TAIR(51»20, _.l_„ IAIP(6 )b32. ...) TAIR(7)b32. t TAIR(8)«32. 

00013 TAIR(9)b20, I IAIR(10)«20. I TAIR(11)b20, i TAIR(13)b10i 

00014 XIN(1)b7,0 I XIN(2)b|,I | XIN(3)>9,3 I XIK(4)b9,8 

. 00015 XIN(5)B9,8_.t.. XINtS)B9.7 X1N(7 )b9,5 | X1NC8)b9,0 __ 

00016 XIN(9)a8.8 I XIN(10)«8,0 I XIN(ll)a7,3 I XIN(12)«7.2 

■ 00017 DO 1 JB1,12 

.00018 i XIN(J)bX1N(J1/24,_._ .... 

00019 NPlTii(i,99} 

00020 99 rORMAT(T30f 'PPKJI bTOTAL ANNUAL COST/DEN(J) Jal, 2< 3' f />T20i 'NHEREl 

00021 2’i/» T)0,*PEN(naSUMMAHON QELEC',/iT)0,'DENta)*SUMHATION OUSE'f/, 

00022 3l30»'DENC3)aSUHHATION OHORK'.///) 

00023 Hib,0i 5 I SIGHAaS>6697E>ll I EPSa,04 I TRa29a| ETARaEPr/100, 

.00024 DO 2 N»li8 _ _ _ _ 

00025 WPITE(),995) XK(K) 

00026 999 rOPHAT(/f IX. *(MDOT*CP/AREA ABSORBER)a',rt0,-4(/> 

. 00027 NPITEHi 100) 

00028 100 rOPHATtn. 'CONC'Tll»'aB»'T21»'PPl(l)'T36,'Pri(2)*T90,*PFI{S)'I66,*DE 

00029 2»DEN(i)»T8l,'DEN(2)'T96f 'DENn)'»/iTll,'aSSSSB*T2lt'8/KWH*T36,'8/KHH 

.000)0 )H'I51( *J/KWH'T66, 'KWH/M2*YR*T81, 'KWH/H2-YR*T96,»KMH/H2-YR',//)_ _ 

00031 IFCICODE ,E0, 2) GO TO 3 

000)2 TCBIF4J, I BETAS, 0024 I TCHAXa441. 1 CCELLaCELLl 

00033 CO TO 4 _ 

000)4 3 ICaTF+5, I BETAa,004l | TCHAXa268. f CCELLaCELL2 

000)5 4 CONTINUE 

00036 TCXsTC + 27), I TFKaTF*273, |TRKaIR + 273, _ 

00037 TYbXK(K) 

00036 XFCTY ,LE, ,1) IDELaS 

00039 IFCTY .CT, ,29) IDEtal0. _ _ _ 

00040 IFCTY ,GE, ,6) IDELb2S 

00041 IFCTY ,GE, ,8) IDELaSB 

.00042 _ CTHS0, _ . _ 

0004) 18 CTHaCTH^IDEL 

00044 . DO 6 JB1,12 

00045 TAIPKsTAIRCJ)+27), 

00046 ICCS0 

00047 19 ICCSICC41 

. 00048 _ IFCICC ,CT. 50) 8TOP_*CONVEBOENCE_PROBIEM8* __ _ 

00049 TCKlBTCK«FCTCX>XINCJ).IAIRK,XKCK))/FPCTCKtXIN(J)iTAXRKfXK(K)) 

00050 IFCABSCTCKl-TCK) ,IT. 1,0E*3) CO TO 1*6 

_00051 .. TCKbICRI 

00052 GO TO 15 

0005) i6 TCbTCK-273, 

_00054 PPal,-BETA*CTC-IP) 

00055 IFCRP ,LE, 0,) GO TO 9. 

00056 GCOOLaXKCK)«CTC>TF) 



rEurcr iN8iB2,rop 


FORTRAN V. 4(218) /KI 


1I>AUC>TS 


19132 PACE 1>1 


000S7 CE{J)»ETAl»XIN(J)«CTH*ETAR«(l,-BETA*(TC-TB)l 

0005t CHORK(J)*,S»U,«TFK/TCK}*QCOOL +QE(J) 

00059 . 8USE(J)»OE{J)+aCOCL _ 

000A0 DO 7 J«l>13 

000it GE(J)sOE(J)*24,*DATH(J] 

00053 OU5E(J5»OUSEf J)»34,«DAyH(J) 

0005) 7 CWORK{J)»OHORK(J)»34,«DATH(J) 

00054 DENi«0, I DEN2a0, | DEN]a0. 

00055 _ DO I J»1M2 

00056 DENIeCENI^OCCJ) 

00057 CEN2«DEN2*QUSECJ) 

00058 J OFNlaOEN)+QWORK(J) 

00059 DEhl«OENl/CTH | o'eN 2»DEN2/CIH j DEN3«DEN)/CtH 

00070 CCOST«CONfETAt«l0, * CCELb/CIH 

0007J A»(l,+XI)«»N' 

00072 ■■(A*XI/(A>l,))aXH 

0007) XNUH«B*CCOST 

tZil* _ _ Prit«XNUH/DENl I PFI2»XNUH/DEN2 ) PFDaXNUH/DEN J 

00075 TI»Pri2+PEI) 

00075 tiRITE().200) CTH,TT,PFIl«PFI2,PFI),DENWDEN3iDEH) 

00077 200 rORMAT(tX,r7,3,3X,F9,3,3Xt6(SX.tPE14.S)} 

00071 S CONTINUE 

00079 irc RR ,GT, 0,) 00 TO 18 

00080 _ 2 CONTINUE 

00081 RETURN 

00083 END 


SUBPROCRAMf CAUED . 

fP 

ABS, r 

BCAIARs'aND ARRATS t NO EXPLICIT DEFINITION • «%" NOT REFERENCED 3 


.F0007 

! 

.F000) 

2 

TAIR 

3 

XK 

?2 

■"#N 

33 

♦ BETA 

24 

• TRK 

30 

♦DENI 

31 

♦ IR 

33 

XIN 

36 

•ICKt 

53 

OUSE 

53 

,F0002 

71 

♦CON 

72 

♦ ICC 

73 

•CEL12 

77 

♦OCOOL 

100 

♦ J 

101 


*A lJ^ ,8000) 121 _ 180002 122 

.80000 (25 aCELM 127 QWORK 1)0 

•XK 142 aTCKAX 16) aCTH 164 

.F0004 )70 ,F0000 171_ *TAIRK 172 

• SIGRA ' 174 ■ '•CCEll 177 ' ♦TT 300 

temporaries 

(PCRtS ,00000 ))0 

PERFCR r MO ERRORS DETECTED ) 


•EFF 

17 

♦DEN2 

20 

♦ IDEL 

21 

♦CCOST 

35 

♦K 

26 

♦ hi 

27 

♦ B 

33 

♦ICODE 

34 

♦ PFI) 

38 

♦ TFK 

57 

♦TF 

70 

%.F0005 


♦PF12 

74 

♦XI 

75 

♦XNUH 

75“ 

♦ PFIl 

103 

QE 

103 

.80004 

117 

♦ ETAR 

123 

.80001 

124 

♦KIAl 

_12S_ 

,F000S 

144 

.F0001 

145 

DATM 

145 

♦TCK 

15S 

♦TI 

155 

♦IC 

157 

♦RR 

173 

♦EPS 

174 

♦DEN) 

-L7»_ 


00001 ))1 


00002 132 


00003 3)3 



SDATt LN8LB2iroft 


roPTKAN V, 4(310) /KI 


ll-AUC-79 


19133 PACE 1 


08001 SUBROUTINE SDATA(BETA>TeHAX,Ht»IN,ETAl,Tr>EFFfTAIR>XK) 

80003 REAL IN 

0000J DINENSION XKt») 

80004 8ICMA«5,669TE>|1 I EP9«,04 I EIARnEFF/100, > TR*3S, 

0000S DO 2 N*ti8 

_00004 NR1TE(3,995) XKtKl _ 

0000T 999 F0RMAT(/,1X>'(MD0T»CP/AREA ABSORBER)a*«F10.4i /) 

00001 NRZTE(3,100) 

. 00009 „._.i80. _ rOPMATtn,'CONC'tll.'TEHP*T2«»'OElEC'I41i'«COOL'T96r'ONORK'T7li*OLt'I84» . 

00010 3T'I0«, 'OUSE'TICl* 'OEtEC/OUSE*Tll6, 'W0UT/0U8E'#/) 

00011 TC>Tr*5, I TAIRK«TAlRf371, 

_0001J TCK*TC*27). _ , 

0001) 4 CONTINUE 

80014 IF( IC ,CE. TCMAX) GO TO 9 

.00015 A»3,»H1»(TC»TAIR)+3.»EP5»SICHA*C(TCK#»4)*(IA1FK**4))+XK(K)»(TC»TFJ 

80016 B«ETA1*IN«(1 .•ETAR«ETAR*BETA«(TC-TR)) 

00017 CTH«A/B I OCOOL»XK(K)«(TC-TF) 

_00016 0E»ETAl*lN«CIH*EIAF*Cli-BEIA*(TC-TR)) _ _ 

00019 NOUT»,5»(1,-(TF+37I,)/TCK)*OCOOL 

80020 OUSEaQE^QCOOL 

_ 80021 OWOPK*«OUI+OE .. . . ... . _. _ 

00022 OLT*QCOOL«WOUT 

8002) FRA1«QE/QUSE I FRA3*NOUT/QUSE 

_00024 _«PlTEt),3001.CIHiIC.QEiQCOOL,awORX.QHiaiUItIPAl.FR.A2 

80025 388 FOPKAT(tX,F7.3>3X,F6.2r6X,7(lX,lPE14,9)) 

00026 TOTC + 20. 

_ 80027 ICK«IC+27),.. . 

80036 CO TO 4 

80039 5 CONTINUE 

000)0 3 CONTINUE 

800)1 RETURN 

800)3 END 


.8UBPP0GRAM8 CALLED 


__SCAlARa ANP.ARRAT8 .t J!»"_».0 EXPL1CIT_DEFIMITI0N._» J^%*_.NOT REFERENCED J 


•TAIR 

1 

• EFF 

3 

XK 

3 

•BETA 

4 

• K 

9 

•HI 

9 

•TR 

.7 

• B 

.18 

•QU8E 

11 

•TF 

13 

•FRA3 

1) 

•QCOOL 

14 

•OE 

IS 

•A 

16 

• NOUT 

17 . 

• ETAR 

30 

•EIAl 

31 

•FRAJ 

23 

<80000 

3) 

•QHORK 

24 

IN 

39 

•TCMAX 

36 

• 9LT 

37 

•eiH 

)0 

. «TCK 

31 

*TC 

33 

«TA1RK 

3) 

•EPS 

34 

• SICMA 

39 




TEMPCRARIES 
isbAlS 104 

_fDA.IA __ t NO JRB0R8 DETECTED J 



MAIN, LNSLB2,rOR 


rORTRAN V. 4(210) /KI 


18»AU6>79 


19123 PACE 1 


00001 

00002 

0000 ) 

00004 

00005 

00006 

00007 

00008 

00009 

00010 
00011 
00012 
0001 ) 

00014 

00015 

00016 

00017 

00018 

00019 

00020 
00021 
00022 

00023 

00024 

00025 

00026 

00027 

00028 
000)9 
00030 
000)1 
000)2 

00033 

00034 
000)5 
000)6 


C- MAIN program for SOLAR CONCENTRATOR CHARACTERISTICS 
DIMENSION XK(8) 

COMMON ErF,CON,CELLl,CELL2*N,XI,XH,ETAl,XK,IF 
REAL IN 

XK(l)ii,01 t XK(2)*.0) I XKM)s, 0 S , XK(4)m,1 I XK(S}«,29 

XX(6)»,5 I XK(7)».75 t XK(8)»,9 

BETA*, 0024 I TCHAX*44l, | Hla,01S | IN*$,0E*i ) ETA1«,T 

TFs60, I Err»15, I TAIR»25, 

NPITE(3,233) 

23) rCRHAT('l'»lX,*GA-AS SOLAR CELLS',//) 

CALL SDATA(BETA,TCMAX,Hl,IN,ETAl,Tr,EFF,TAIR,XK) 

__ 1>RIIE(),2)4) _ 

J34 rORMATC'l', IX, 'SILICON SOLAR CELLS',//) 

BETA*. 0041 I TCMAX*268, 

CALL SDATA(BETA,TCMAX,Hl,IK,ETAl,TF,Err,TAIR,XK) _ . _ 

NPIIE(3,23)) 

EFF*18, I BETA*. 0024 J TCHAX*441, 

CALL SDATA(FETA,ICMAX,Hl,IN,ETAl,Tr,ErF,IAIR,Xl(} 

ErF»15, 

CON*100, 


C 



CELLl*i00000a 
N*20 I XI*. 08 
CALL BLOCXl 
CALL 8LOCF2 
VARY CELL COST 
CELLl*80000, t 
CALL BLOCXl 
call BLOCX2 
CELLl*20000, > 
CALL BLOCXl _ 
CALL BLOCX2 
CELL1*8000, I 
CALL BLOCXl 
CALL BLOCX2 
STOP ' NORMAL 
END 


CELL2*S000, 
I XM«,03 

CELL2*5000. 
CELL2*20000, 
CELL2*)000, 
END' 


common blocks 

/.COMM./(*2n 

EFF _ *0 _ CON ♦! CELLl +2 CELL2 _ ♦) N +4 

XI is ■ XH" ♦S' EtAl +7 ■ XK ♦IB It +20 

subprograms called 

SDATA BLOCXl BLOCK2 

SCALARS ■anD''APRAYS~ t **« Ny EXPLlCIT“bE>WlTION i%«'“H 0I RErERENCES ') 

•TAIR L _ *BEIA 2_ IHl... _IN _ 4 fTCMAX 5 
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